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ABSTRACT 


This study presents the development and methodology^ for 
development of full -authority implicit model- following and 
explicit model -following optimal controllers for use on 

helicopters operating in the Nap-of-the Earth (NOE) environment. 
The controllers were designed based on NOE handling qualities 
requirements and the longitudinal dynamics of the conventional 
AH-1G and the AH-1G fitted with a simplified longitudinal 
- auxiliary propulsion system. Pole placement, input -output 

frequency response, and step input response were used to evaluate 
handling qualities performance. The pilot was equipped with 
velocity -command inputs. 

A mathematical/computational "trajectory optimization" method 
was employed to evaluate the ability of each controller to fly NOE 
maneuvers. The method determines the optimal swashplate and 
thruster input histories from the helicopter's dynamics and the 
prescribed geometry and desired flying qualities of the maneuver. 
Minimization of a cost function which incorporates these dynamic, 
geometric, and flying qualities constraints resulted in the 
optimal control histories. From the control law for each 
controller, these optimal helicopter control histories were used 
to — "back out" the pilot input histories. The helicopter state 
histories and pilot input histories were used to evaluate the 
controllers in terms of their abilities to meet the flying 
qualities criteria and to minimize pilot workload; pilot workload 
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was measured by the integral cost of stick input rates. The 
method avoids the time and cost constraints, and biases., 
associated with piloted simulation studies. 

Three maneuvers were investigated for both the implicit and 
explicit controllers with and without auxiliary propulsion 
installed: pop-up/dash/descent , bob-up at 40 knots, and 
glideslope. -The pop-up and bob-up clearly showed the advantages 
that the use of auxiliary propulsion has in reducing the time to 
complete a task and in providing increased longitudinal 
acceleration/deceleration while not requiring excessive pitch 
attitudes . 

The explicit controller proved to be superior to the implicit 
controller in performance and ease of design. The design syt hcsis 
showed that, in most cases, it was necessary only to weight the 
diagonal elements of the state and control weighting matrices; 
guidelines for selecting these elements is included in the study. 
Pilot inputs using the explicit controller were smooth, decoupled, 
..and. in proportion to the desired horizontal and vertical 
velocities of the helicopter. 
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Chapter I 


INTRODUCTION 


1.1 Background and Motivation for this Study 

Today's battlefield helicopters need to be highly maneuverable. 
In particular, the maneuvers that must be performed for successful 
Nap -of- the -Earth (NOE) tasks require a high level of helicopter 
agility, good handling qualities, and reasonable levels of pilot 
workload. These agility and handling qualities criteria often 
dictate both the required response characterisitics of the control 
system, and the required acceleration potential of the helicopter. 

Such specifications can often only begin to be realized through a 
full -authority controller, one that can utilize the full control 
power of the helicopter. Additional requirements to rapidly 
transition from one point to the next suggests the need for higher 
acceleration/deceleration levels, which may make the inclusion of 
auxiliary propulsion mandatory. 

Motivated by these considerations, the aim of this study was to: 

1. define the ''operational effectiveness" criteria for NOE 
operations ; 

2. investigate suitable controller structures for 
full -authority control of the helicopter which meet the 
established criteria and which are sufficiently easy to 
design; 

3. investigate the advantages incurred by using 
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auxiliary longitudinal propulsion in NOE operations; 

4. introduce a mathematical/computational method for 
evaluating the performance of a controller and auxiliary 
propulsion system for NOE applications which does not employ 
human pilots or pilot models. This method is called the 
trajectory optimization method. 

1.2 Definitions Of Handling Qualities and Agility 

At the onset, it is necessary to define and distinquish several 
terms which are used throughout this work. 

Handling qualities are "those qualities or characteristics of 
an aircraft that govern the ease and precision with which a pilot 
is able to perform the tasks in support of the aircraft role.” 
[14] Agility is "a measure of the ease with which a helicopter 
can change its state.” [10] 

There is an important distinction between handling qualities 
and agility. Handling qualities are allied with the dynamic 
response of the helicopter, and as such, handling qualities 
criteria define the helicopter's response to a given control 
input; quantitatively, they define how the controller should 
perform in terms of damping ratio, bandwidth, control 
sensitivities, pole/zero locations, control force gradients, et 
cetera. Agility is more closely allied with translational motion. 
Quantitatively, agility can be measured in terms of particular 
aspects of the helicopter's state time histories, most notably its 
accelerations. Handling qualities criteria, agility, and pilot 
workload are the three measures used to evaluate the performance 
of the helicopter and controller over a prescribed mission. These 
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three measures collectively constitute what may be called the 
mission's "operational effectiveness." 

1 . 3 Overview of this Study 

This study addresses the application of full -authority optimal 
controllers and auxiliary propulsion to helicopter NOE operations. 
Chapter II describes the NOE environment and gives both a 
qualitative and quantitative account of the desired flying 
qualities and handling qualities within this environment. In view 
of these criteria, the advantages of using auxiliary propulsion 
are discussed, as is the preferred controller command structure. 

Chapter III presents the formal development of the explicit 
model -following and implicit model -following controller 
structures. They are developed using the linearized equations of 
motion of the helicopter. Chapter IV explains how the controller 
and auxiliary propulsion system are evaluated for a given 
maneuver. A trajectory optimization method is used to compute the 
optimal helicopter input control histories. The pilot input 
histories are found from these control histories and the control 
law; pilot workload is measured from the pilot input histories. 
Chapter V presents a detailed description of the mathematics and 
computational solution for the trajectory optimization method. 

The implicit and explicit model -following controllers are 
applied to the AH- 1G helicopter. The selection of the controller 
gains for both of these controller structures, at the hover trim 
condition, is presented in Chapter VI. The examples in this 
chapter are given for both the conventional helicopter and for the 
helicopter equipped with auxiliary propulsion. 
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Three NOE maneuvers are investigated in this study: 
bob-up/dash/descent, pop-up at 40 knots, and glideslope. The 
optimal control input histories, computed using the trajectory 
optimization program, are shown in Chapter VII along with the 
histories of the state variables. Finally, in Chapter VIII, the 
pilot input histories are shown for each of the three NOE 
maneuvers, using both of the controller structures, and for the 
helicopter with and without auxiliary propulsion. The 
"operational effectiveness" of all of these helicopter/controller 
configurations is evaluated in this chapter. 



Chapter II 


THE NOE ENVIRONMENT 

2.1 General Attributes of NOE Missions/Maneuvers 

The roles of the pilot and helicopter vary in the battlefield 
environment. But in general, NOE maneuvers are employed chiefly 
to avoid detection by the use of cover. [10,14] The helicopter 
moves discretely from one point to the next. Maneuvers may be 
elusive, as during air-to-air combat; they may be operational, as 
in transporting troops and/or cargo; or they may be of the 
scout/attack type, for example air-to-air combat and anti-tank 
missions. The latter usually employ target search and weapons 
delivery; many employ navigation. A generalized list of necessary 
attributes for flying NOE maneuvers is: [10,12,13] 

* high dash speed 

* good longitudinal/vertical acceleration/deceleration 

* small turn radius 

* good handling qualities and controllability 

* low pilot workload 

2.2 Typical NOE Maneuvers 

The literature commonly defines a list of 15-20 high and low 
speed maneuvers descriptive of NOE flight: [10,12,13,15,19,28] 
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Low speed 

* precision hover 

* taxi 

* lateral jink 

* rearward flight 

* bob -up 

* dash 

* slalom 


High speed 

* box pattern (four turns) 

* longitudinal accel/decel 

* pul 1- up /push- over 

* turning approach 

* straight approach 

* lateral accel/decel 


* quick-stop 

* liftoff and landing (including ship deck liftoff/landing) 

* pop-up (zero forward velocity) 

* hover about a point/masked hover 


This work is concerned with the longitudinal degrees of freedom of 
the helicopter. We therefore investigate tasks that include the 
bob -up and pop-up, dash, quick- stop, longitudinal 
acceleration/deceleration, and straight approach maneuvers. 
Actually, our purpose in selecting this group of maneuvers is 
twofold. This group reflects representative longitudinal 
maneuvers, and it also reflects those maneuvers for which 
auxiliary longitudinal propulsion is perceived to be most useful. 
Sample trajectories are computed and shown in Chapter VII. 

The bob-up is usually associated with a quick ascent from 
behind tree cover or other terrain cover; pop-ups are like bob-ups 
but begin from hover. The dash is a sprint maneuver typically 
from fifty feet to a half mile. The dash is most commonly 
employed to rapidly transition from one point to another; as such, 
it is often accompanied by a quick- stop and descent to hover. 
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Other longitudinal maneuvers require high speed flight over longer 
distances. A typical bob-up maneuver at high speed is obstacle 
avoidance, for example in "hurdling" a bridge while flying low 
over a river. 

2.3 Desired NOE Flying Qualities 

Two flying qualities issues stand out as being very important 
in NOE flight. These are the pitch attitude constraint and 
acceleration/deceleration potential . 

2.3.1 Pitch Attitude Constraint 

It is important in NOE maneuvers to maintain a desirable pitch 
attitude. Most often this attitude is close to zero. These small 
pitch attitudes are very important for carrying out secondary 
operations, armament and navigation in particular: [10,13] 

O For weapons delivery, maintaining near-constant pitch 
attitude is crucial for radar "locked on" tolerances. 

O Large nose-up and nose -down pitch attitudes are undesirable 
for navigation and target search. 

O Large pitch forward attitudes raise the the tail rotor 
above tree -line, making it easier to be seen and heard. 

O Pitching must be small for the pilot to maintain visual 
cues . 

O There are dangers of tail rotor ground strike and autorota- 
tion with excessive nose-up attitudes. 

Though there is not universal agreement on exactly what are the 
permissible attitudes for NOE maneuvers, the literature suggests 
that pilots prefer nose-down attitudes in the 0° to 15° range and 
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nose-up attitudes from 0° to 10°. [10,12,19] 


2.3.2 Acceleration/Deceleration Requirements 

Time is the most important agility parameter. Most tasks are 
performed in minimum time, or, at the highest possible 
accelerations and decelerations. But because conventional 
helicopters require large attitude changes for acceleration and 
deceleration, the acceleration/deceleration levels are restrained 
because of the pitch attitude constraints outlined above. In a 
survey of Army, Navy, and Air Force pilots [10,19], the ability to 
accelerate and decelerate -- at higher rates and with less pitch 
attitude -- was sighted as the issue in which they would like to 
see the biggest improvement. 

Figure 2-1 plots acceleration versus airspeed for a typical 
helicopter. The plot shows that acceleration potential decreases 
with airspeed. Particularly for battlefield operations, pilots 
prefer to have higher acceleration potential across the airspeed 
spectrum. The current ABC [25,29] and Tilt-Rotor [28] aircraft 
increase the acceleration potential at higher airspeeds. These 
programs, and the studies conducted in References 10, 16, 24, 25, 
and 29 define the desired potential acceleration level, across all 
airspeeds, at 0.3 to 0.5g. 
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Figure 2-1 Available acceleration vs airspeed 

(Figure adapted from References 24 and 25) 

2.4 Potential Benefits of Auxiliary Propulsion 

In light of these flying qualities criteria, we can predict the 
improvement in agility that would result using an auxiliary 
longitudinal thruster on a helicopter operating in this 
environment. The addition of an additional X-force via a 
longitudinal thruster would increase available acceleration 
levels, so that the acceleration potential of the helicopter looks 
like the upper line in Figure 2-1. And, because this X-force is 
thruster- generated rather than rotor-generated, large pitch 
attitudes are not required to achieve greater 
acceleration/deceleration levels, so that pitch attitude 

constraints should be met. These predictions are substantiated in 
Chapters VII and VIII. 

2.5 Handling Qualities and Controller Selection 

As one might expect, there is a direct connection between the 
type of maneuver the pilot is expected to perform and the type of 
command input which the pilot prefers to have. For the 


9 


longitudinal degrees of freedom, the flight (i.e. state) variables 
are the helicopter horizontal and vertical velocities, and pitch 
attitude and pitch rate. Numerous studies have investigated 
longitudinal controllers with various command inputs. The two 
basic configurations are velocity command configurations and 
attitude/attitude-rate command configurations. From the 
literature available, it is determined that for the types of NOE 
maneuvers investigated in this study, pilots prefer velocity 
command controllers. The following paragraphs summarize some of 
the work found in the literature and explain the rational for 
choosing the velocity command controller configuration over 
alternative configurations . 

The pilot's control over the full longitudinal and 
lateral-directional motions of the helicopter typically consists 
of four command inputs; pitch, roll, yaw, and vertical command 
inputs. For longitudinal motion, these inputs reduce to the pitch 
and vertical command inputs. Although more (or fewer) pilot 
inputs may be defined, two inputs is the most efficient for 
control within this horizontal-vertical plane. Ideally, one input 
should control vertical motion, and one input should control 
horizontal motion. [6,11,12,13,15,28] 

Most NOE tasks require the pilot to govern the speed and 
position of the helicopter; other tasks, notably air-to-air 
combat, require attitude control. It is, of course, logical to 
equip the pilot with command inputs that will best enable him to 
perform the tasks at hand. For the maneuvers of interest in this 
study, which are velocity/position based, we find that pilots 
prefer a rate (velocity) command controller, with possible 
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modifications. The results from a literature review of controller 
structures for NOE tasks, outlined below, show that velocity 
command systems are best for maneuvering, and attitude command 
systems are best for precision hover tasks. 

References 17, 28, and 31 conclude that pilots prefer velocity 
command over attitude command (AC) for tracking tasks. Typically, 
both velocity command and attitude command systems have a vertical 
velocity command as one input; for the other input, the velocity 
command system uses a horizontal velocity command and the attitude 
command uses a pitch attitude command. This preference also holds 
for low speed precision pointing tasks; reference 19 found that an 
AC controller degraded the helicopter's handling qualities during 
precision pointing tasks. The study conducted in that reference 
also found that the addition of attitude hold (AH) to the velocity 
command controller results in a high level of agility, as measured 
by achievable rates and accelerations. 

For chase and aquisition, reference 31 finds that pilots prefer 
the velocity controller. For these types of tasks, the large 
stick inputs necessary with the AC controller are found to be 
objectionable; the AC controller is preferred for tasks where 
precision maneuvering is critical. This reference, and references 
11, 13, and 15 note that the best handling qualities are achieved 
when the two velocity commands in the the velocity system are 
decoupled; decoupled contols is highly desirable for good handling 
qualities. Hence, the study conducted in this thesis uses an 
velocity controller. 

With the control structure thus defined, it was then possible 
to set some quantitative response criteria for the velocity 
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controller that would result in good handling qualities. Many of 
the defining handling qualities criteria used in this work were 
taken from the Army's Advanced Digital/Optical Control System 
(ADOCS) program. [11,15] This control system is being designed to 
provide satisfactory handling qualities for attack helicopters. 
The ADOCS report itself includes an extensive literature search to 
define applicable handling qualities criteria for attack 
helicopters . 

At the onset, the two velocity command inputs should be 
decoupled, and should produce responses that closely approximate 
the outputs of two first-order systems. Specifications for the 
two first order time-constants vary strongly according to pilot 
and task, but the literature generally agrees that the accepted 
range for the vertical response time constant is from 1 to 4 
seconds, while the horizontal response time-constant ranges from 2 
to 5 seconds. [11,13] When a first-order response cannot be 
achieved, a second- order response should have a damping ratio 
equal to or greater than 0.7; the minimum acceptable damping ratio 
is 0.5. 

The input frequency band of concern does not vary greatly among 
the references [11,13,15,18,20], and thus the band was taken to be 
from 0 to 1 Hz (6.26 radians/second). In addition, the transfer 
function phase response should follow that of a first-order 
system, with a suitable time-constant as specified above. Too 
much phase shift can result in pilot induced oscillations (PIO) . 
[ 20 ] 

For the velocity- command controller, the vertical velocity 
command input sensitivity should range from 13 to 16 ft/sec per 



inch of stick, for the vertical response time-constants listed 

above. [11,13]; the horizontal velocity input sensitivity ranges 

from 12 to 16 ft/sec per inch, for the related time - constants . 

This range is again due to the variation in tasks from which these 

data were taken. These inputs should not produce pitch 

2 

accelerations in the helicopter which exceed 0.69 rad/s . 

The handling qualities criteria are summarized in Table 2-1 
below: 



horizontal 

velocity 

command 

vertical 

velocity 

command 

first-order 

time-constant 

2-5 sec 

1-5 sec 

damping ratio 

>0.7 

>0.7 

minimum 
damping ratio 

0.5 

0.5 

bandwidth 

0-6 rad/s 

0-6 rad/s 

control 

sensitivity 

12-16 ft/s/inch 

13-16 ft/s/inch 

maximum 
pitch rate 
sensitivity 

.69 rad/s/inch 

.69 rad/s/inch 


Table 2-1 

Handling qualities criteria for velocity- command controller 
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Chapter III 


CONTROLLER DEVELOPMENT 


3.1 Linearized Equations of Motion 

Aircraft performance analyses and controller development are 
commonly accomplished using the linearized equations of motion of 
the vehicle about a trim point. These linearized equations can 
very closely approximate the total nonlinear equations of motion, 
and are much easier to use in application. In addition, 
linearized numerical data for many aircraft is readily available. 
The controller development methods and trajectory optimization 
procedure of this study utilized the helicopter's linearized 
equations of motion. All numerical data was derived for the AH-1G 
helicopter, which typifies helicopters that are used for the roles 
discussed in this study. 

Appendix A reviews the procedure and assumptions used in the 
linearization of the helicopter's nonlinear equations of motion. 
It is shown below that the helicopter's total state variables are 
the sums of the trim states and the linearized (perturbation) 
states. Appendix B contains the numerical (linearized) stability 
and control derivatives for the AH-1G at twenty- two trim points 
from reference 1. A realistic representation of the helicopter 
states at all points and times in space requires that these 
derivatives be scheduled between trim points. Appendix B 
discusses how these data are scheduled according to flight 
condition. 
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The speed, 

position, 

acceleration, 

and 

orientation of 

the 

helicopter are 

described 

in terms of 

its 

nominal (trim) 

and 


perturbation state variables. Figure 3-1 shows the total and 
perturbation states of the helicopter. Because of the small angle 
assumption used in the linearization of the equations of motion, 
the total state variables U, W, Q, and 0 are the sums of their 
trim values (U , W , Q , 9 ) and perturbation values (u, w, q, 6 ). 

o o o o 

Since perturbation positions (x and y) are measured along the 

trimmed body axis, their contribution to the total inertial 

referenced positions must take into account the trim attitude 9 . 

r o 

Therefore , 



Figure 3-1 Longitudinal Body Axis System 


3.2 Design Objectives 

The objective of this study was to design a controller for 
helicopters that perform NOE tasks, or other demanding tasks, 
utilizing the full available control power of the helicopter. 
Such a controller is a "full -authority" controller, one which is 
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capable of exploiting the full control potential of the 
helicopter, rather than one which limits the controls to less than 
their full range of travel (as is often done in 
stability-augmentation systems) . There are numerous controller 
structures which can potentially fit this bill. For controller 
design using the helicopter's linearized equations of motion, 
several popular and promising controller structures are those 
referred to as linear -quadratic optimal controllers^ 

A suitable controller structure is one that satisfies the 
prescribed input/output response criteria (handling qualities 
criteria) . But it is also of concern that the controller be 
sufficiently easy to design , and that it is adaptable to changing 
flight conditions (i.e., easy to schedule within the flight 
envelope). Linear- quadratic implicit and explicit model- following 
optimal controllers prove to meet these requirements. 

Pilots have shown favorable responses in studies conducted on 
aircraft using model -following controller structures. Reference 
18 shows that pilots rated model-following better than response 
feedback mode control, producing respective Cooper-Harper ratings 
of 3 versus 3.5. The pilots described a sense that the 
model -follower controller was "locking on" when it was engaged. 
Reference 17 also notes pilots' preferences for model- following 
over response feedback control because of the ability of the 
model -following system to provide inter-axis’ control to compensate 
for undesirable cross-coupling affects in the vertical and 
horizontal axes. 

The ADOCS program [11,13] used a model -following concept to 
achieve the desired command responses. In that study, classical 
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control methods were used to design feedback paths that met the 
stability requirements. However, based on the desired 
command/response characteristics, a command reponse model was used 
to provide the desired response. It achieved good pilot ratings, 
falling within Level 1 of the Cooper-Harper ratings scale. This 
feedforward command augmentation and shaping performs the same 
function as the feedforward gains of the model- following 
controllers presented below. The feedforward structure provides 
control mixing (to achieve decoupled modes) and prefiltering, but 
does not affect the level of feedback stabilization. The 
model- follower structure also facilitates easy flight condition 
scheduling . 

3.3 Linear-Quadratic Optimal Controllers 

The method employed in the design of optimal controllers is to 
minimize a predetermined cost function, subject to the system's 
dynamic constraint, for the purpose of generating an optimal 
control law, i.e., an optimal relationship between the states and 
controls, that gives the resulting closed-loop system the desired 
control/response characteristics. The general method is to 
define, from the cost function and system dynamics, the equations 
that establish necessary and sufficient conditions for minimum 
cost. [2,4] The necessary conditions are derived by equating the 
cost function's linear sensitivity to zero* through a procedure 
known as the Calculus of Variations. The cost function is 
typically comprised of a final state penalty and an integral 
penalty function of the state and control. Parameters of the 
integrand, or Lagrangian (£) , define the nature of the optimizing 
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solution. The cost function is written as: 


J - <£[x(t f ) ,t £ ] + J £ £[x(t) ,u(t) ,t]dt 


(3-1) 


The system dynamics: x - f [x(t) ,u(t) , t J 

and initial condition: x(t ) « x 

o o 


(3-2) 

(3-3) 


must be satisfied while J is being minimized. To assure adherence 
to this "dynamic constraint," these dynamics are adjoined to the 
cost function via an adjoint (or influence) vector l(t). 


J - <£[x(t f ) ,t f ] + 


£{■ 


2[x(t) ,u(t) ,t] 


+ 1 (t) {f [x(t) ,u(t) ,t] -x(t) Udt 




(3-4) 


f [x(t) ,u(t) , t] -x(t) is adjoined to the integrand because it must 
be satisfied over the entire interval <t ,t>. And because it 

o’ f 

equals zero when satisfied, the numerical value of J is not 
changed. l(t), then, acts as an influence function, and expresses 
the cost function's sensitivity to dynamic effects. 

The Hamiltonian, H, of the integrand is defined as: 

J([x(t) ,u(t) ,t] - Jf[x(t) ,u(t) ,t] + l T (t)f [x(t) ,u(t) ,t] (3-5) 

and J can be rewritten as 

J - rf[x(t £ ),t f ] + J t £ |3f[x(t),u(t),t] - l T (t)3c|dt (3-6) 

(Explicit reference to time is dropped here; variables at the 
final time t f are denoted by the subscript . 

Using the Calculus of Variations, three necessary conditions are 
derived from the cost function. These are: 



(3-7a) 


1 . 

2 . 

3. 


0 - (^) T 
K du ; 


i - (- ^) T 

' 3x' 


s - <S-> T 


(3-7b) 

(3-7c) 


f 

These equations are also known as the Euler -Lagrange equations. 
When satisfied, the Hamiltonian is stable for infinitesimal 
control variations. These conditions are therefore local, rather 
than global criteria. The sufficient condition is that 3f >0. 
For linear -quadratic cost functions (defined below) , H >0 is 
guaranteed to be satisfied. Therefore, the formal mathematics for 
this sufficient condition are not pursued here. 

Inspection of Equations 3-2 and 3-3 and Equations 3 -7b and 3 -7c 
shows that this optimization problem is two -point boundary value 
problem; the state integration constants are specified at t^ and 
the adjoint integration constants are specified at t f . 

For linear-quadratic controllers, the dynamics are linear and 
the Lagrangian consists of weighted quadratic norms. A 
quadratic, time -invariant, final state penalty takes the form 
^[x ] - X fPf X f* and t * ie Euler -Lagrange equation becomes l f * 
P f x f . However, the integration interval is taken from 0 to t f , as 
t f approaches infinity. This means that the final state penalty 
is presumably insignificant. [4] With this in mind, and recalling 
that 1 and x are adjoint and therefore Hold over the entire 
interval <0,«>, we drop the final time subscript for the third 
Euler -Lagrange equation, and now have a more suitable, yet wholely 
correct, condition that 1 - Px. 

The sufficient condition, ^ uu >0, is guaranteed to be satisfied 
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for quadratic cost functions as long as the weighting matrices are 


positive-definite. The elements of the weighting matrices are 
parameters that are chosen to meet the closed- loop response 
criteria; a unique optimal controller results for each uet of 
weighting matrices chosen for the giver system dynamics. The 
relationship of weighting matrices parameters and control 
objectives is not direct nor is parameter selection intuitive. 
However, linear-quadratic model-following optimal controllers 
prove to be easy to develop because one usually needs only to 
adjust the diagonal elements of the weighting matrices, and 
perhaps an off-diagonal element for obviously coupled variables. 


3.4 Implicit Model -Following Controller 

The system (in our case, the helicopter) is 


dynamics : 


x - Ax + Bu 


where for our application 


u 

w 

q 


and 


u 


Is 


described by the 
(3-8) 

(3-9) 


and A and B are the stability and control derivative matrices, 

respectively. Note that we include an auxilary longitudinal 

thrust control, T (inclusion of T does not alter the derivation). 

The model, whose state response (x ) to control inputs (d) we 

ro 

want the system (helicopter) to emulate, is expressed as 


x - Fx + Gd (3-10) 

m 

where x is the plant state vector as given above and d is the 
vector of (pilot) inputs. For the velocity-command controller, d 
contains horizontal and vertical commands: 



d 


u 

com 

W 

com 


(3-11) 


This is an implicit model -follower structure because the error 
-- in this case, the rate error -- between the system and the 
model is penalyzed IruuxLe the cost function, rather than as part 
of the control system in the closed- loop system dynamics. [9,21] 


The cost function is: 

J « — J j\x-x ) T Q(x-x ) + u T Ru + 1 T (Ax+Bu) - l T xj-dt (3-12) 
where x and x are nxl, Qisnxn, uismxl, and R is m x m. 

m 

Substituting for x and x from Equations 3-7 and 3-8, 


ii:{ 


(Ax+Bu-Fx-Gd) Q(Ax+Bu-Fx-Gd) + 


u T Ruj-dt 


-J |(x T [A-F] T +uV-d T G T )Q([A-F]x+Bu-Gd) + u T Ruj-dt 




x T (A-F) T Q(A-F)x + u T (B T QB+R)u + d T G T QGd 


T/1 t-> \ - rt 


+ 2x (A-F) QBu - 2x (A-F) QGd - 2d G QBuj-dt (3-13) 


Using the following definitions, 


W - (A-F) T Q(A-F) W -= (A-F) T QB W - g t qb 

xx XU du 


W - (A-F) QG 

xd 


W - BQB+R 


W - GQG 

dd 


(3-14) 


the cost function becomes 


1 00 f 

J -= — f ix T W x + u T W u + d T W d + 2 x T W u - 2x W d 

2J 1 xx uu dd xu xa 

- 2d T W u + 1 T (Ax+Bu) - l T xidt 
du J 


(3-15) 
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and the Hamiltonian is 


n - -U T W x + u T W u + d T W d + 2x T W u - 2x T W d 

2 1 xx uu dd xu xd 

- 2d T W u + 1 T (Ax+Bu)1 (3-16) 

du J 

Employing the three optimality conditions, we have 

1. 0 - (|^) T - W u + W T x - W T d + B T 1 (3-17) 

OU UU XU du 

, 

2. i « (^) T - -U x - W u + W d - A T 1 (3-18) 

OX XX XU xd 

At this point we observe that the integrand contains the 
variables 1, u, x, and d. Therefore, we revise the algebraic 
expression for 1 to include effects of the model input d. This is 
done by adding a linear inhomogeneous term to the (linear) 
homogeneous expression for 1 from the previous section [22]: 

3. 1 - Px - Sd (3-19) 

P is n x n and S is n x r. Solving Equation 3-15 for u yields 

U - -W’ 1 (W T x - W T d + B T 1) (3-20) 

uu xu du 

Using Equation 3-19 we substitute for 1 in Equation 3-20 to obtain 
the optimal control law in terms of the desired variables x and d: 


- -w ' 1 

(W T x - W T d 

+ B T Px - 

B T Sd) 



uu 

xu d u 




- -w " 1 

T T 

(W +B P)x + 

W _1 (W T + 

B T S)d 

(3- 

-21) 

u u 

X u 

u u du 




- -C x 

+ C d 



(3- 

•22) 


The expressions for P and S are found from the dynamic equation 
for 1. Sustituting u from Equation 3-21 and 1 from Equation 3-19 



into Equation 3-18, we have 

i - -W X - W [-W' X (W T +B T P)x+W" 1 (W T +B T S)d] + W d -A T (Px-Sd) 

xx xu uu xu uu du xd 

- f-W +w w' x w T +W w' 1 B T P]x 

XX xu uu xu xu uu 

+ f-W W' X W T -W W" X B T S+W +A T S]d (3-23) 

xuuudu xuuu xd 

Taking the derivative of Equation 3-19 with respect to time, and 

noting that the best estimate of d is zero (since it is 
unpredictable [22]), we have 

1 - Px + Px - Sd (3-24) 

Substituting for x from Equation 3-8 and then for u from Equation 
3-21, and collecting like terms, Equation 3-24 becomes 

i - [P+PA-PBW* X W T -PBW _ 1 B T P]x 

u u XU u u 

+ [PBW'V +PBW" 1 B I S-S]d (3-25) 

u u d u u u 

Equating the right hand sides of Equations 3-23 and 3-25, and 
recognizing that, for minimum cost, the resulting equation will 
hold true irrespective of the values of x and d, we generate two 
separate equations by equating the coefficients of x and the 
coefficients of d in these equations: 

P - -P[A-BW'V ] - [A-BW'V ] T P 

UU X U U U X u 

+ [-W +W W‘ X W T ] + P[BW* X B T ]P (3-26) 

XX xu uu xu uu 

S - -[A-BW'V ] T S + P[BW'V]S 

uu xu u u 

+ P[BW'V ] - [W -W W'V ] (3-27) 

uu du xd xu uu du 

Equation 3-26 takes the form P - -PM -M T P + PNP + K and is 
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recognized to be a matrix Riccati equation in P. The steady- state 
P matrix is solved according to the diagonalization method [26] 
outlined in Appendix C. Once P is determined, the ^teady-state 
solution for S (S - 0) is determined from Equation 3-27 using 
simple matrix algebra. 

Equation 3-22 is the control law for the implicit 
model -follower. Substituting this expression for u into the 
open- loop dynamics, the closed- loop dynamics are 

x-Ax + Bu*=Ax + B(-C x+C d) * (A-BC )x + BC d (3-28) 

1 2 v V 2 

Figure 3-2 shows a block diagram of the open- loop helicopter, and 
the closed- loop helicopter with feedback gains C on the measured 
states x, and feedforward gains from the pilot inputs d. 


closed-loop helicopter 



Figure 3-2 Implicit Model Controller 


3.5 Explicit Model -following Controller 

The open- loop dynamics of course do not depend on the 
controller structure, and are again given by Equation 3-8. The 
idealized model dynamics are 
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(3-29) 


x - Fx + Gd 

m m 


where x is the model state vector and d is the pilot control 

m 

input vector. 

The linear -quadratic cost function for the explicit 
model -follower penalizes control excursions and the error between 
the actual state and the model state. [2,21] 


(x-x ) T Q(x-x ) + u T Ru + l T (Ax+Bu) - 1 xldt (3-30) 

m m J 

[x T Qx - 2x T Qx + x T Qx ] + 1 T (Ax+Bu) - l T x)dt (3-31) 

1 m m m J 

Because there are two distinct states, we could define an 
augmented state which includes both x and x^. The cost function 
using the augmented state is 




J 



X + 


T 

uRu 


+ 1 (Ax+Bu) - 


i T 

1 x 


:jdt 


(3-32) 


where 

X - 

X 

and 

i 

Q - 

Q -Q 



X 

L m j 



,-Q Q . 


Minimization of the cost functions given by Equations 3-31 and 
3-32 yield dimensionally different optimality equations. The 
resulting control law in each case, however, yields identical gain 
coefficients for the state and command vectors x and d. 
Derivation of the control law using Equation 3-32 is exactly 
analagous to the derivation for linear quadratic regulators (LQR) ; 
this method is not presented here. Instead, the derivation using 
the "unaugmented " cost function, Equation 3-31, is used as 
perhaps it best shows how the control law depends on X m as well as 
x and d. 
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The Hamiltonian is 


H « -Hx T Qx - 2x T Qx + x T Qx + u T Ru + 1 T (Ax+Bu) 

2 I m m m 


The first two optimality conditions are 


0 - ^> T - Ru + fiTl 


• 8H T T 

2. 1 - (- g) 1 - -Qx + Qx -A~l 

OX m 


(3-33) 


(3-34) 

(3-35) 


As we did in the implicit model -follower development, we assume a 
linear relationship among 1 and the variables of the integrand 
(other than u) : 

3. 1-Px + Px -Sd (3-36) 

12m N 


From Equations 3-34 and 3-36, the control law is 


u - -R'Vl - -R' 1 B T (P x + P x - Sd) 

1 2 m 

= -R _1 B T P x - R" 1 B T P x + R' 1 B T Sd (3-37) 

12m V 7 

*= -C x - C x + C d (3-38) 

1 2 m 3 ' 7 


To find expressions for P^ , P and S, solve for 1. 
Substituting Equation 3-36 into Equation 3-35: 

i *= -Qx + Qx - A T (P x + P x - Sd) 

m 1 2 m 

- [-Q - A T P ]x + [Q - A T P ] x + [A T S]d ■ (3-39) 

1 2 m 

Taking the derivative of Equation 3-36 with respect to time (and 
recognizing that the best estimate of d is zero) , we have 

i-Px+Px+Px + P x -Sd (3-40) 

1 1 2 m 2 m . 
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Substituting Equation 3-8 for x and Equation 3-29 for x , and then 

m 

substituting Equation 3-37 for u and collecting terms, we have 
i - [P A+P +p br" 1 b t p ] x + [P br" x b t p +P F+P ]x 

L lll 1 J L 1 222m 

+ [P BR' 1 B T S+P G-S ] d (3-41) 

1 1 2 J 

Equating the coefficients of x, x , and d from the right hand 
sides of Equations 3-39 and 3-41 yields the expressions for P^, 

P , and S : 

2 


P -= -P A 

- a t p + p br" x b t p 

- Q 

(3-42) 

l i 

11 1 



P «= -P F 

- a t p + P br" x b t p 

- Q 

(3-43) 

2 2 

2 1 2 



S = P G - 

T -IT 

AS + P BR B S 


(3-44) 


2 1 


Equation 3-42 is a matrix Riccati equation; the steady-state 

P^matrix is solved according to the diagonalization method [26] 

described in Appendix C. Once P^ is determined, steady- state P^ 

is solved from Equation 3-43 via a Kronecker product method [30], 

since P^ premultiplies and postmultiplies other (constant) 

matrices. With P and P known, steady-state S is solved 

12 

algebraically: 

S - [A T - P BR* X B T ] " X P G (3-45) 

The block diagram for a helicopter equipped with an explicit 
model -following control law (Figure 3-3) shows that the ideal 
model acts as a prefilter of pilot inputs d to the open- loop 
helicopter . 
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Figure 3-3 Explicit Model Controller 


3.6 ^ther Controller Structures 

Other linear-quadratic optimal controller structures were 
investigated. In particular, the proportional- filter and 
proportional- integral structures were considered. [2] However, 
choosing suitable weighting matrices parameters for these two 
structures proved to be a very laborious assignment. Because of 
this, and because the response characteristics of these 
controllers showed no improvement over the model -following 
controllers, they were eliminated as possible controller 


structures . 









3.7 Gain Scheduling 

The controller gains are derived for the linearized dynamics of 
the helicopter about a trim point. As such, they are only 
applicable for flight within a region close to this trim 
condition. For the controller to be effective throughout the 
entire flight envelope, it became necessary to derive gains for 
trimmed flight conditions throughout the helicopter's entire 
flight envelope, and to schedule the gains between these discrete 
trim points. Scheduling is typically done by expressing the gains 
as a function of flight condition, e.g., airspeed, advance ratio, 
and dynamic pressure. [8] Reference 7 suggests a means to select 
the flight variables and the gains which require scheduling by the 
use of correlation analysis. Any scheduling method should be 
chosen to provide accurate gain values throughout the entire 
envelope, while not being overburdening in terms of computation 
time and data, or cost. 

The data available for this study limits the number of trim 
data points to twenty- two. This means that twenty- two sets of 
gains have been derived for each controller structure investigated 
(see Appendix D) . But this limited data set makes it feasible 
only to schedule the gains in terms of the most obviously 
influential variables, the helicopter's forward and vertical 
speeds, U and W. And it is impossible to heuristically argue what 
functional dependence the gains have with U and W between trim 
points (unlike what can be done for the stability and control 
matrices -- see Appendix B) . Therefore, the gains are linearly 
interpolated and extrapolated using the trim data set. 

Figure B-l shows the spread of the data (trim) points in the 
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U-W plane. In this study, a data set contains the value of every 
gain element at each of these twenty- two points for each 
controller. For each [U,W] pair sampled throughout a trajectory, 
the gain elements are scheduled by first interpolating with 
respect to U, and then with respect to W. For flight conditions 
that put [U,W] outside of these data points, the data is simply 
extrapolated out to the [U,W] point. 



Chapter IV 

CONTROLLER AND AUXILIARY PROPULSOR EVALUATION 

4.1 Method of Evaluation 

In Chapter I, we expressed the operational effectiveness of the 
helicopter/control system as a measure of its ability to perform 
its prescribed mission. In the context of this study, we are 
specifically concerned about the abilities of the controller and 
auxiliary propulsor to meet this goal. System performance, or 
operational effectiveness , is described by three primary groups 
of criteria: handling qualities and flying qualities, both of 

which are outlined in Chapter II for NOE operations, and pilot 
workload. The handling qualities criteria are evaluated 

coincident with the development of the controllers (Chapter VI) , 
with one exception, controller sensitivity, discussed below. 
Flying qualities and pilot workload must, of course, be evaluated 
over a specific maneuver. This study used a "trajectory 
optimization" method for evaluating these two categories of 
criteria. Use of this method avoids the traditionally subjective 
rating method for evaluating pilot workload and contoller 
sufficiency. 

4.1.1 Trajectory Optimization Method 

The trajectory optimization method computes, for a maneuver of 
specified time, geometry, and flying qualities, the optimal 
helicopter control histories, i.e., the motion of the helicopter's 
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physical controls (swashplate and thruster input) throughout the 
trajectory. The process is purely mathematical, and is performed 
computationally. Flying qualities and the geometry of the 
maneuver are incorporated by specifying the helicopter's state 
vector at various times throughout the trajectory. Excursions 
from these ’'constraints" are penalyzed according to their 
importance, and, using the dynamic equations of motion of the 
helicopter, the best (optimal) control input history is 
determined. A pilot implicitly performs exactly the same process 
in attempting to fly a specified maneuver. [10,12,17,19,33] 
Chapter V describes the trajectory optimization process and 
algorithm in detail. 

4.1.2 Flying Qualities and Pilot Workload 

Using the optimal control histories derived from the 
optimization, the flying qualities and pilot workload evaluations 
were then performed as follows: 

1. ) If an optimal control history could be found that forced 
the helicopter to fly the prescibed trajectory with the prescibed 
flying qualities, then we concluded that the helicopter was 
capable of achieving these qualities. An inability to meet the 
criteria (e.g., minimum transition time) suggested a deficiency in 
the helicopter (e.g., lack of sufficient acceleration potential). 

2. ) From the control law for each controller, we "backed out" 
the pilot input time histories. These histories enabled us to 
evaluate the controller for satisfactory control sensitivity 
(criteria tabulated in Section 2.5) and decoupled response, and 
enabled us to compare the pilot's "workload" over that trajectory. 



The "workload" was a quantitative measure derived from the pilot 
input-rate histories; the process used here is explained fully in 
Section 4.3. 

The premise for using this method is straightforward, and the 
method's advantages are readily seen: 

1. ) This method does not employ human pilots or a simulator. 
It is purely computational, but much simpler compared to a 
simulator, though equally effective for preliminary analyses of 
controllers and the effects of auxiliary propulsion. In addition, 
the method certainly does not require the cost and time 
commitments of simulation studies, nor does it depend on the 
availability of pilots with comparable training. 

2. ) Taking the pilot out of the loop permits separate and 
unbiased evaluation of the helicopter attributes (notably agility 
via available control power) and the control system attributes. 
The inherent ability to satisfy the geometric and flying qualities 
criteria throughout a trajectory, with acx^eptofa&e pilot workload, 
is essentially a function of the features of the helicopter, and 
not of the pilot. But, because human pilots are very adaptable 
[14], an ill-equipped helicopter may be made to fly as desired at 
the expense of the pilot's workload. Consequently, the pilot 
rating system makes it impossible to separate out pilot workload 
associated with "inadequacies" of the helicopter, from workload 
associated purely with the functioning of the control system. 

In the trajectory optimization method, the trajectory's 
geometry and flying qualities are easily quantified, and are 
exactly the same criteria which the pilot attempts to meet when 
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performing the maneuver. But by specifying the physical control 
input histories, and not pilot input histories, we make the 
helicopter fly as close to these constraints as possible, and we 
can therefore evaluate the helicopter's attributes independently 
of the pilot and the control system. 

3. ) Pilot workload is effectively described using a simple 
quantitative measure. We quantify pilot workload by observing 
that workload is directly related to the amount which the pilot 
must exercise his controls. That is, workload is directly related 
to the "area under the curve" of the pilot input-rate histories 
(we use the integral of the norm of the input rates -- see Section 
4.3). 

Reference [10] suggests that a quantifiable measure of pilot 
workload can be related to the amount of time that the 
helicopter's pitch attitude exceeds specified upper and lower 
bounds. This may be a good measure, but in the method of this 
study, we force the helicopter to remain within these bounds. 
Therefore, this type of workload, which is really associated with 
the control power of the helicopter, can be disregarded, so that 
our valuation of pilot workload seems to give a better appraisal 
of the controller itself. 

4. ) The control histories extracted from the trajectory 
optimization program are for the open- loop helicopter. 
Determination of these optimal control histories is therefore 
independent of the controller. The inputs must come from the 
control system, but are completely independent of the control 
system <Jsuictu/ie . To evaluate the performance of any control 
system structure, we take the optimal input histories, and the 
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control law for that controller, and back out what the pilot 
input histories to the controller must be to generate the 
controller output histories, which are the optimal Input histories 
to the helicopter. This means that, for each maneuver, the 
optimal histories need to be computed only once, and then any 
number of controller structures can be evaluated. 


4.2 Pilot Input Histories 


4.2.1 Implicit Model Controller 

Extracting the pilot input history, d(t), for this controller 
involves nothing more than simple matrix algebra. From Chapter 
II, the control law is u(t) *= -Cx(t) + C d(t) . For a given 
trajectory, u(t) is the (optimal) history of swashplate inputs and 
thruster inputs (if a thruster is included), x(t) is the state 
history, i.e., the trajectory that was prescribed and which 
results from the input history, u(t) , and C and C^ are the gain 
matrices. The control law equation is solved for d(t) : 


d(t) - C 2 [u(t) + C^x( t) ] 


(4-1) 


d(t) is the history of pilot control input commands. This study 
assigns the pilot two controls, vertical and horizontal velocity 


command: 


w 


u 


Therefore, when the thruster is employed, 


there are three helicopter controls, and C^ is not square. The 

solution for d(t) is overdetermined, and we must use the left 

pseudo - inverse of C : C L (C T C ) 1 C . C and C are scheduled 

r 2222212 

according to the total flight velocities, U and W. 
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It is also necessary to scale the command vector, d, so that 
the pilot inputs are within an acceptable range. The obvious 
choice is to scale W to the range of the collective stick, and 

com 

U to the longitudinal cyclic stick range. These scale factors 

com 

in no way affect the stability or response characteristics of the 
control system. The factors form a 2 x 2 diagonal matrix, 


D 


8 / w 0 

c, range com, range 


which simply scales d by 


P / u 

! Is, range com, max I 

premultiplying the feedforward matrix C (the left pseudo- inverse 
is now a weighted left pseudo- inverse) ; see Figure 3-2. 

For the thrust augmented helicopter, the helicopter has three 
control degrees of freedom while the pilot has only two. Using 
the pseudo- inverse to solve for the pilot histories effectively 
was a "minimum norm" fit of d(t). Therefore, it was necessary to 
check that the state history that resulted from pilot inputs to 
the closed-loop system (x - [A-BC^x + BC^d) matched the state 
history determined in the trajectory optimization process. For 
the trajectories and controllers investigated in this study, the 
state variables of both histories (i.e., the integrated 
closed- loop history and optimized history) matched within one 
percent of each other over the trajectory interval. This was 
expected because of the good low frequency input -output responses 
of the controllers. Had these two histories not been consistent, 
a suitable weight for a weighted left pseudo- inverse 
transformation that resulted in matching state histories would 


have been chosen. 



4.2.2 Explicit Model Controller 

The prefilter in this controller makes the solution for d(t) 
slightly more complicated than the solution for the implicit model 
controller. It is necessary to simultaneously integrate the model 
dynamics while calculating the pilot input history. 

The control law is u(t) - -C x(t) - Cx(t) + C d(t) . Solving 

1 2 m c 

for d(t) : d(t) - C L [u(t) + C x(t) + C x (t) ] 

3 1 2m (4-2) 

However, in this equation x (t) is unknown. It must be calculated 

m 

from its dynamic equation, 

x(t) = Ax(t) + Gd(t) (4-3) 

m m 

But since d(t) is present in this dynamic equation, Equation 4-2 
must be solved while Equation 4-3 is being integrated. The 
procedure is as follows: 

The state and model initial conditions are known: x(t ) - x(t ) = 

o mo 

x . u(t) and x(t) are known for all t. Therefore, we can find 

o 

d(t ) from Equation 4-2. Then x at the next step, i.e. x(t +1) , 

o m mo 

is found by integrating Equation 4-3. Then solve Equation 4-2 for 

d(t +1) , and so on up to and including time t . In general, the 
o f 

procedure is 

1. Solve Equation 4-2 for d(t ) using the initial conditions. 

o 

2. For n discrete time entries in the time history, do the 
following for k * 1 to n-1: 

integrate Equation 4-3 to get x(t ) 

m k+1 

solve Equation 4-2 to get d(t ) 

k+1 

Because a fourth-order Runge-Kutta method is used for 
integration, step 2 of the above iterative procedure is actually 
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done four times for each integration time step. And, because the 
gains are dependent on flight condition, they are calculated 
whenever step 2 is performed. 

When d is overdetermined , i.e. when C has more nonzero rows 

3 

than columns because there are more physical controls u than pilot 
commands d, the left pseudo- inverse of C 3 is used in Equation 4-2. 
And, the pilot command vector, d, is scaled to the conventional 
collective and cyclic stick ranges. 

The same discussion of the left pseudo- inverse and state 
history matching of Section 4.2.1 applies to the explicit 
controller . 

4.3 Assessment of Pilot Workload 

As stated in Section 4-1, pilot workload is assessed in terms 
of control sensitivity, decoupled response to commands, and the 
pilot command inp 1 ' rates. The first two are easily understood, 
and used without explanation in Chapter VIII. 

Intuitively, the more the pilot has to move the sticks, and the 
higher the rates of movement, the more work he is doing to fly the 
aircraft. The amount of stick movement, then, is a good 
quantitative measure of pilot workload. The amount of work 
performed in moving the sticks is the integral of the stick input 
rate -- actually, to avoid assigning a negative work value for 
negative rates, we used the norm of each (scalar) rate history. 
Thus, this measure of workload is the "cost” of the rate input, 
using unit weighting: 

r t 

Pilot workload ~ [ f ||d(t)||dt (4-4) 



This measure of workload depends on the time and nature of the 
maneuver. Therefore, it best serves as a means of comparing 
controllers for the oame maneuver. 

This measure was also applied to the helicopter control rate 
histories (collective, cyclic, and thruster rate histories). 
These work "costs" are used as normalizing measures for the 
controller inputs , and also to compare U vs B workloads and 
V vs ^ workloads. 

com c 
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Chapter V 


THE TRAJECTORY OPTIMIZATION METHOD 

This chapter is comprised of three sections. The first 
presents a general description of the formulation and minimization 
of the cost function used in the trajectory optimization process. 
In the second part, this general trajectory optimization method is 
applied to the helicopter problem to determine the optimal control 
histories for helicopter trajectories; individual terms of the 
cost function are physically interpreted in the context of the 
helicopter trajectory optimization problem. These two sections 
are mathematically formal; optimized sample trajectories are given 
in Chapter VII. The third part of this chapter describes the 
algorithm used for the optimization. 

5.1 Formal Mathematics 

We begin with the dynamics of the plant and initial conditions 
of the state: 

x - f[x(t),u(t),t] (5-1) 

x(t ) - X (5-2) 

O O 

where x and u (for our purposes) are the helicopter state and 
control vectors, respectively (see Chapter III). 

Initially, we define a cost function in the manner used in 
designing optimal controllers. [4] The cost function J penalizes 
state, state-control, and control excursions throughout the 



interval. The integrand of "penalties" is the Lagrangian, £. 

J - f f £[x(t) ,u(t) ,t]dt (5-3) 

t 

O 

Since we seek the control history u(t) which yields a stationary 
minimium of this cost function, subject to Equations 5-1 and 5-2, 
we adjoin the dynamic equation to the cost function, thus 
including it in the minimization. Note that when the dynamic 
equation is satisfied, this additional term in the cost function 
equals zero, thus the cost function is numerically unchanged. 

J = J f ji?[x(t) ,u(t) ,t] + l T {f [x(t) ,u(t) ,t] - x(t))jdt (5-4) 

O 

In many applications, one often desires to specify state 
conditions at the end of the interval and at times along the 
interval. Mathematically, this is accomplished by including 
interior point penalty functions, 0 , and a final point function, 
<f > . Incorporating these penalties (also termed constraints), the 
redefined cost function becomes 

i 

J « *[x £ (t),t f ] + Z^[x.(t),t.] 

+ J t f |£[x(t),u(t),t] + l T {f [x(t) ,u(t) ,t] - x(t)}jdt (5-5) 

where i is the number of interior point constraints. 

Finally, it is useful (and realistically imperative) to have 
saturation constraints on the controls. These control constraints 
are expressed by the function 

C[u(t),t] < 0 (5-6) 

As with the dynamic constraint, the control constraint is included 
in the minimization of the cost function by adjoining it to the 
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integrand of the cost function; the adjoint vector for control 
constraint is m(t) . The redefined cost function becomes 

J - <j> [ x ( t ) , t ] + 20. [x. (t) ,t. ] + f f {if[x(t) ,u(t) ,t] 

I I 11 1 Jt 

O V 

+ l T {f[x(t),u(t),t] - ift)} + m ft)C[u(t),t]|dt (5-7) 

m(t) is an influence vector, and each element of m takes its value 
according to whether or not its corresponding control is 
saturated : 

r > 0 C - 0 

m (t) \ 1 (5-8) 

1 ^ - 0 C < 0 

i 

This means that the product of each control/adjoint element pair 
is zero, or mC - 0, so that the cost function is numerically 
unchanged. The cost function now incorporates all the desired 
penalties and constraints. 

A Calculus of Variations procedure is used to derive the 
conditions for optimality. [4] The objective is to find those 
conditions which will cause the cost function to be at a 

stationary minimum. Therefore, we begin by taking the first 
variation (the differential) of J and equating it to zero. 

Explicit reference of the variables to time is dropped here. And, 
for simplicity in notation, only one interior point constraint is 
used; the results are easily generalized to include multiple 
interior point constraints. 

It is convenient at this point to define the Hamiltonian, H, of 


the cost function: 



The Hamiltonian has the same form as the one used for optimal 
controllers, with the addition of the control constraint term. 


Then J 


[x ,t ] + 0 [x ,t ] 

n f f J 1 i i 



l T x + m T C)dt 


(5-10) 


dJ — (f^dt + f^dx) + £ dt + (|^-dt + 4— dx) 

''at f ax 't-t t-t £ at i ax t-t 

£ f i 


+ £ dt + r f {(|^ + i T §| + roT |^) 53 

t-t i Jt 1 dx 3x OX 

i o v 

,b£ -Tflf , Tac, f -T.O. 

+ ^ + 1 an + m dt )Su - 1 5 x j dt 


(5-11) 


A few remarks: This equation contains no terms for the 

differential changes in f and x because differential changes in f 
exactly cancel with differential changes in x due to their 
equality, Equation 5-1. The differential recognizes variations in 
the state, the control, and the derivative of the state -- there 
are no control derivative terms *• at a fixed time, t. 

f T ♦ 

We now simplify the integral term I 1 fixdt. For any 

o 

variables y and v, Jydv - yv - Jvdy. Applying this to 

fl T Sxdt over the interval <t , t >, and taking into 

** of 

account discontinuities in x at the interior point, this gives 




t- 


J t f l T Sxdt - l T 5x| t i - J t A lfixdfc 

i o o 

+ l T fix| f - f 1 i T 5xdt 
■t + Jt + 


(5-12) 

But since x is continuous in the interval, the right hand side of 


Equation 5-12 becomes 


l T 6x| t i + l T fix| t f + - J t f + i T «xdt - 1 T 

o i i 

+ l T (t f )6x(t f ) - l T (t i+ )6x(t i+ ) 


( ti -)5x(t.-) 



f • T 

1 fixdt 


l T (t )fix(t ) 

O O 

(5-13) 


We eliminate fix terms by solving for fix from the following 
equations for the differential, dx: 


dx(t £ ) -= 
dx(t £ ) - 


6x(t f ) + x(t f )dt £ 
r fix(t i -)+x(t.-)dt. 
fix(t +)+x(t +)dt 


(5-14a) 
(5 -14b) 
(5 -14c) 


or, fix ( t f ) - dx(t £ ) - x(t f )dt £ 

6x(t -) — dx(t ) - x(t -)dt 
i i i i 

6x(t + ) *= dx(t ) - x(t +)dt 

i i i i 


( 5 - 15a) 
(5-15b) 
(5-15c) 


Using Equations 5-15, and noting that the variation in the initial 
conditions is zero, the right hand side of Equation 5-13 becomes 


[l T (t.-)-l T (t.+)]dx(t. ) - l T (t.-)x(t.-) + l T (t +)x(t. + ) 

ill ii i i 

t 

+ l T (t f )dx(t f ) - l T (t f )x(t f )dt f - J t f i T fixdt (5-16) 

O 

Finally, substituting for x using the dynamic relationships 


x(t-) - f , x(t +) -=f , and x(t ) — f , the workable form 

i t- i t+ f t 

i i £ 


j: 


f T * 

for | 1 fixdt becomes 

t 

O 


I f l T fixdt ~ [l T (t -) -l T (t +) ]dx(t ) - l T (t j -)f i dt 
J t i i i * 


i t - i 
i 


+ l T (t +)f dt + l T (t )dx(t ) - l T (t )f dt 
i t + i f f f t f 

i f 


- 


fixdt 


(5-17) 



4 

Substituting this result into the cost function differential, 
collecting like terms, and writing all suffixes as subscripts, 


dJ 


- J(|£ + £ + l T f )dt + - 1 T )dx) 

1 at t f t f t f y f ax t £ J t-t f 

+ J (M + + 1 T f - 1 T f )dt + (|^-l T +1 T )dx|- 

Vat t £ v t.- v v ' i ax v v JW 

*j::{ 


.sr.Tsr iSc^iT., . ... iac, ,.1 . 18 , 

< S ^‘ 1 a ? 1 >‘* + ‘ST 1 Si™ aS ) 5 u j dt (5 18> 


To achieve J - 0 (minimization of the cost function) , we choose 
1 such that the coefficients of variations in the state, control, 
and time -- dx, fix, 6u, and dt -- equal zero over the entire 
interval <t ,t >. Doing this yields the following conditions for 
optimality: 


1 . 

1 t 8<f> 

t “ ax 

f 



(5-19) 

2. 

•t d£ .rdf Tac a« Tac 

1 ”'ax' 1 ax' m ax“‘ax ax 



(5-20) 

3. 

” su au au au au 



(5-21) 

4. 

0-1^ + 2 + l T f 

at t-t l t Vf 

f 



(5-22) 

5a . 

i - i 

t - t + ax 

i i 



(5-23) 

5b. 

0-|* + £ +l’f - 1 T £ ' 

8t ,.t l , V V V V 

i 



(5-24) 

The sufficient condition -- that the second 

variation 

of J to 

change 

s in x, u, and t must be positive -- is 

not 

treated here. 

It is 

noted below that this condition holds 

for 

the 

helicopter 


application problem. 
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Equation 5-19, together with Equations 5-20 and 5-21 without 
the m term, are the classic Euler -Lagrange equations for fixed 
end-time optimal control problems. Equations 5-23 and 5-24 are 
the optimizing conditions when the end time and interior point 
times are allowed to vary; i.e., these are the optimizing 
conditions for open end- and interior- time problems. Equation 
5-23 defines the discontinuity in the adjoint history l(t) at the 
interior point time t . Equation 5-19 specifies a final time 
adjoint condition, and because x is specified at t , this is a two 

o 

point boundary value problem and necessitates an iterative 
solution (Section 5.3). 

5.2 Interpretation for Helicopter Trajectories 

The cost function used for the helicopter trajectory 

optimization is linear -quadratic -- the helicopter dynamics are 
linear and the cost function penalties are quadratic norms of the 
state and control weighted by conformable positive -definite 
weighting matrices. Positive-def initiveness of the weighting 
matrices assures that the stationary point (found via cost 
function minimization) is a minimum stationary point. 
Mathematically, pos it ive - definitiveness assures that the second 
variation of the cost function is positive at the stationary 
point. The method of using quadratic weighting matrices with 
linear dynamics is analagous to the method ‘used in developing a 
linear-quadratic regulator. 

This section discusses each term of the cost function in the 
context of the helicopter problem. 



5.2.1 Dynamic Constraint 

The dynamics of the helicopter can be viewed as a dynamic 
constraint which restricts the state accelerations to a given 
function of the state and control. The dynamic equation is x - Ax 
+ Bu. For the purpose of trajectory optimization, these dynamics 
are augmented in three ways: (1.) to include perturbation 
positions; (2.) to include the addition of auxilary propulsion to 
the helicopter; and (3.) to account for control rate limits. 
These three items are discussed separately. 

5. 2. 1.1 Augmenting the dynamics to include change in position 
Just as the perturbation rates and attitude, [u,w,q,0], must be 

added to their respective trim values, [U W Q 9 ] , so too must 
the perturbation positions, [x y] , be calculated and added to 
their trim values, [X Y ] . Perturbation variables u, w, q, 6 are 

O O 

defined in a body axis system (Figure 3-1), so that perturbation 
positions from trim must be resolved to the trim axis: x - ucosfl 
and y - wcos0. The dynamics are augmented to include these 
perturbation states: 


u 

w 

q 

6 

x 

y 




9 in the cos0 terms makes these augmented dynamics "slightly" 
nonlinear . 

5. 2. 1.2 Inclusion of Auxilary Propulsion 

In the current longitudinal framework, the two physical 
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controls of the helicopter are the swashplate controls: 
collective and longitudinal cyclic. Addition of an auxilary 
propulsive device would institute another control, which is 
included in the dynamics as follows: 


u 
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(5-26) 

It is assumed that the auxiliary propulsion system provides purely 
longitudinal thrust. Thus, as Equation 5-26 shows, the only 
nonzero control derivative in the third column of the augmented B 
matrix is the derivative which directly affects horizontal 
acceleration, u. This is obviously a very simplified approach. A 
much better understanding of the auxiliary propulsion unit and its 
aerodynamic and inertial effects on the helicopter is needed 
before its effects on any of the stability and control derivatives 
can be quantified. As for X^, there is little data available on 
stability derivatives for longitudinal thrusters on helicopters or 
VTOL aircraft. A value for is best determined by sizing it 
according to the desired acceleration level. This is discussed in 
Sections 2.3.1 and 6.1. 

5. 2. 1.3 Control Rate Limits 

As noted in Section 5.1, the cost function constrains the 
controls to fall within their saturation limits. This is, of 
course, a realistic necessity since a helicopter (or any other 
system) has physical limits on its control power. But it is also 



necessary to constrain the control oatea, since the hydraulic 
actuators impose a delay on the command input. These limits are 
included in the optimization as follows. 

Consider the scalar control u, and assume it behaves as a 
first-order system: 

u - -— u + — u (5-27) 

T T com 


For a unit step input, 


i 


the response is 


u - (1 - 


-t/r . 

e ) 


(5-28) 


and 



so u 

max 


_1 

r 


(5-29) 


This means that the control is rate -limited by ±— . Employing this 

T 

method of rate- limiting has two particular advantages. Firstly, 
this first-order lag realistically models the actuator dynamics of 
the swashplate (Figure 5-1). The actuator lag time-constant is r 
in the equations above. Secondly, rate -limiting the controls in 
the dynamics means that a rate -limit constraint does not have to 
be explicitly included in the cost function. Including such a 
constraint would add more complexity to the cost function and 
increase the optimization computation time. If the constraint was 
explicitly included, it would be appended to the cost function 
with an influence vector in a fashion similar to that done for the 
control constraint. [4] This would, however, necessitate 
back- differencing the control history to compute the control rate. 
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I 



physical 



other controller 


inputs ( e . g . , 


feedback 


and prefilter signals) 


Figure 5-1 Actuator lags 


Each of the three controls is rate -limited by a first-order 
system with an appropriate time - constant . These control dynamics 
are appended to the dynamics of Equation 5-26. The physical 
controls, $ , B , and T now become part of the state vector and 

c Is 

the commanded controls now serve as the control vector. 
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We now let x - Ax + Bu represent this augmented system. This is 
precisely the dynamic constraint, x - f[x,u], used in the cost 
function. 

5.2.2 Flying Qualities Constraints 

The Hamiltonian penalyzes state and control excursions from 
zero throughout the flight. We can loosely say that the 
Hamiltonian quantitatively characterizes the desired "flying 
qualities" of the flight. Consider the quadratic Hamiltonian 

W[x,ul - -x T K x + x T K u + -u T K u (5-31) 

The first term penalizes excursions of the state. Term three 
penalyzes control excursions and term two penalyzes state-control 
excursions. Excursions of a particular variable are penalyzed by 
assigning a positive weight to the corresponding element(s) of the 
appropriate weighting matrix. Because the Hamiltonian is 
integrated, and because the optimization procedure works to 
minimize the Hamiltonian, those state and control elements which 
are weighted must necessarily be driven toward zero. That is, 
their excursions (from zero) are minimized. Motions of state and 
control elements with zero weight are unrestricted -- though of 
course the dynamic constraint and other applicable constraints 
still apply. 

For example, it is often desired to restrict the pitch rate and 
pitch attitude of the helicopter throughout the trajectory. In 
this case, the [3,3] and [*,«] elements of the matrix are 
assigned positive weights. Intuition and experience dictate that 
it is usually only necessary to assign values to the diagonal 


51 


elements. In all cases, however, it must be remembered that the 
weighting matrix of the Hamiltonian -- the one formed by 
augmenting x and u -- has to be symmetric positive-definite. 
Determining the value of a weighting element is a matter of 
"engineering judgement." One develops a "feel" for sizing these 
weights according to their observed influence in the resulting 
optimal trajectory. 

One final point about the weighting matrices. .-'This analysis 
uses constant, time- invariant matrices. Matrices that are a 
function of time are probably not realistically necessary. 
Matrices that are a function of the state and/or control 
perturbation or nominal -- may be effective. However, if a 
weighting matrix is a function of a perturbation variable, the 
optimality conditions that result from the partial derivative of 
the Hamiltonian with respect to this variable will be slightly 
changed. 

At this point, with dynamic constraints and flying qualities 
constraints, the cost function is 

J ~ f f /-x T K x + x T K u + u T K u + l T (Ax+Bu) - l T x 
J t I 2 x xu u 

o v 

5.2.3 Interior and Final Point Constraints 

Most NOE maneuvers are associated with ^achieving a specified 
velocity or position at a specified time, or within minimum time. 
Therefore, it is necessary to be able to specify the state vector 
at a point in the trajectory. This ability equates mathematically 
to constraining the state to equal a specified value at this 
specified time. 


h 


(5-32) 



Again, quadratic constraints are used for the interior and 
final point constraints. Denoting x as the desired state vector 
at the final time t f , the error between the actual state and the 
desired state at t is x -x . Using the weighting matrix K , we 

f f fd 1 
can construct a quadratic penalty of the final state error: 


Similarly, for (multiple) interior points, 


(5-33) 


20 


1 

[x(t ),t ] - -E(x -x ) T K (x -x ) 

2 i id i i id 


(5-34) 


Including these constraints, the cost function is now 

i 

J - -(x -x ) T K (x -x ) + iz(x -x ) T K. (x. -x ) 

2 f fd f f fd 2 i id i l id 

+ f f i-x T K x + x T K u + u T K u + 1 T (Ax+Bu) - l T xldt 
J t I 2 x xu u J 


(5-35) 


5.2.4 Control Constraints 

The swashplate controls and the thruster control have minimum 
and maximum saturation limits. These limits cannot be exceeded, 
and a hard constraint must therefore be imposed on the controls . 
Adjoining the control constraint to the cost function assures that 
the constraint will not be violated. 

The control limits are simply 

u < u < u 
min max 

or, 6 < 6 < 6 

c c c 

min max 

B < B < B (5-36) 

Is Is Is 

min max 
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T < T < T 

min max 


Written functionally as a constraint, 


C[u(t) ,t] 


u-u 


-u+u 


< 0 


(5-37) 


This constraint is adjoined to the cost function with m(t) exactly 
as the dynamic constraint was adjoined with l(t) . Note that the 
dimension of m is twice the dimension of the control, u. 

With the addition of the control constraint, the cost function 
takes its completed form: 


J - -(x -x ) K (x -x ) + -Z(x -x ) K (x -x ) 

2 f fd f f fd 2 i id' i i id 

i 

+ f f /-x T K x + x T K u + u T K u + 1 T (Ax+Bu) - l T x + m T cldt 

Jt x XU u J 

(5-38) 


5.2.5 Optimality Conditions 

Applying the optimality conditions derived in Section 5-1 to 
this linear-quadratic cost function yields the following 
optimality conditions for the helicopter optimal trajectory 
problem : 


1. 


1 - = K (x -x ) 

t 3x f f fd 

f 


(5-39) 


o i / 3Jf T3C.T „ v T, 

2. 1 - (-T- - ® "o — ) — -K x - K u - A 1 

OX OX X XU 


(5-40) 


3. 0 - H “ - K T x + K u + B T 1 + m T C (5-41) 

u O U OU xu u u 


where for three controls, C - diag[i 1 1 -1 -1 - 1 ] 
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4. 0 - 


+ £ + l T f 

at t , t t_ 


t- 1 t 


- (x T K x + x T K u + u T K u) + l*(Ax+Bu) (5-42) 

x XU U t t t 


f f 


d6 


5 “- l , - - l . . ♦h* 1 . ♦ + K i < v i n> 


(5-43) 


5b. o - + 2 + 1 T f - 1 T f 

at t t - t - t + t + 


t-t i 

i 


i i i i 


- (x T K x + x T K u + u T K u) + 1 * (Ax+Bu) _- 1 * £Ax+Bu) 

* XU H t t t t - t i t i + 

(5-44) 


1 


5.3 Optimal Trajectory Computation Algorithm 

The numerical solution method used to solve the trajectory 
optimization problem is the second-order gradient, or steepest 
descent, method. [4,32] This is an iterative method that improves 
the control history, u(t) , on each iteration. Improvement 
connotes minimizing the cost, thereby coming closer to satisfying 
the optimality conditions. 

The following flow diagram shows the logistics of the steepest 
descent algorithm used for the trajectory optimization. A 
discussion of component functions of the algorithm follows the 
flow diagram. 
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Inputs 

trim conditions 
interior and final points 

weighting matrices - interior point 

- final point 

- state 

- control 

- state/control 

time interval and integration increment 
perturbation initial conditions 
initial control history guess 
actuator time -constants 
auxiliary thrust on/off flag 


Control Power 

determine perturbation control limits 


Integrate states x(t) 
integrate forward in time using initial 
guess of control history 


Calculate cost J 
record this cost 


Is cost minimized? 


Integrate adjoint l(t) 
integrate backward in time 


Calculate adjoint m(t) 
for control saturation 
















Figure 5-2 Flow Diagram of Trajectory Optimization Algorithm 
5.3.1 Inputs 

Each interior and final point is specified by the time at which 
it occurs, and by the values of each state element at that time. 
Obviously, t f is the time for the final point. Initially, the 
interior point times must be guessed. Equations 5-42 and 5-44 are 
the optimality conditions which dictate how these initial guesses 
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should be changed for minimum- time problems. That is, these 
conditions are used as a gradient to adjust the interior and final 
point times until the new times result in trajectory histories 
which satisfy the equalities of Equations 5-42 and 5-44. In 
practice, finding the times for minimum- time problems is 

difficult. Rather than employ these optimality conditions in an 
outer -loop iteration around the optimization problem, once an 
optimal (but not minimum time) trajectory has been- computed, the 
"tweaking" of these times is left to the user. The rationale is 
that the outer -loop iteration would have to choose a fraction of 
the interior point time gradient with which to perturb the 
interior point time. Choosing the correct fraction is very 
difficult, and a bad choice could mean considerable time to 
convergence, if convergence occurs at all. The user performs 
essentially the same outer- loop process by observing the completed 
trajectory and control histories, and deciding how much, if any, 
he wants to try to improve the trajectory by changing the interior 
point times. Mechanically, this means respecifying the interior 
point times and rerunning the optimization. The time to rerun the 
optimization is short, since presumably the new interior point 
times are not radically different form the existing times, so that 
the new optimal control history does not differ greatly from the 
existing history. The most efficient method is for the user to 
change the interior point times. 

There is virtually no limit to the number of interior points 

that can be specified -- the limit is the number of time steps 

taken in the interval, minus the ones for t and t , or 

o f 

(t -t,.)/ (integration increment) - 1. In practice, it would never 
o r 


be necessary to specify this many interior points. State elements 
which are oat to be constrained at interior or final points are 
assigned an arbitrary value, but their corresponding weighting 
matrix elements rruuxt be zero. 


The integrand weighting matrices (K^, K^, K^) must be 

conformable with the state and control vectors. They must also be 
selected such that they form a Hamiltonian whose second 
derivative with respect to the control (H ) is positive. In 

uu 

other words , 


K K 

X XU 




must be symmetric positive semi-definite. 


In many applications, K and K are null matrices, so that 
cannot be null. Elements of K are usually small; the identity 

u 

matrix is a good choice for K^. 

Trim conditions are [U (ft/sec), W (ft/sec), Q (deg/sec), 8 

o o o o 

(deg)]. Perturbation initial conditions should always be zero. 

The time interval <t ,t> is from t -0 to t . For minimum- time 

Of O f 

problems, determination of minimum t is best left to the user via 
the method discussed above for interior point times, for the 
reasons given. 

The integrations in this study are performed digitally and are 
therefore discrete. Integrations are performed at 40 Hz, or every 
.025 seconds. The lower limit should be small enough to properly 
calculate the fastest dynamics. In this study, the fastest 

dynamics are the actuator dynamics, which have time constants of 
. 08 seconds . 
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The better the initial estimate of the control history, the 
less time the optimization will take. The example of Section 7.4 
shows how the trajectory converges to the optimal history from the 
initial history. The initial history can be quite dissimilar from 
the optimal history and still converge, as long as the pitch 
attitude is not driven to excessive amplitudes (near 90° and 
larger). And, it of course makes sense to specify control values 
for the initial estimated control history which are within the 
saturation limits. This initial control history mu&t have time 
entries consistent with the specified time interval and 
integration increment. 

The perturbation control limits are determined from the trimmed 
stick positions and the absolute authority limits of each control. 
The perturbation control, in terms of stick movement, is the 
inches of stick available from the trimmed position to these 
minimum and maximum limits. Perturbation stick movements are 
determined by finding the trimmed stick position (in inches) using 
the trimmed swashplate angle (in degrees) and the total stick and 
swashplate travels. For the longitudinal cyclic, this gives: 


r B 

IB 1 


r IB 

B 

Is 

1 S 


1 s 

1 s 

trim 

min 


trim 

min 

B 

B 


IB 

IB 

Is 

Is 


Is 

1 s 

^ max 

min J 

deg 

^ max 

min J 


(5-45) 


inches 


or, solving for IB^ 

trim, inches 


r IB 


IB 


Is Is 

trim min 


Is 


trim, inch 


Is 


1 s 


min J 


x fB - OB 1 + B 

^ Is Is J Is 


deg 


min in min , i n 

(5-46) 


The perturbation stick travel is then 
B - flB - B ] 

Is '•Is Is J 

max, inches max trim inches 

B — flB - B ] 

Is '•Is Is J 

min, inches min trim inches 


(5-47a) 


(5-47b) 


Similarly, for the collective perturbation stick travel: 

e - fe ■ e ) (5-48a) 

c ^ c c J 

max, inches max trim inches 

B C “ C e c - e c ) (5-48b) 

min, inches min trim inches 

Note that the total permissible perturbation stick travel is equal 
to the total stick displacement; the trim position acts merely to 
define the perturbation "zero reference" position. Because the 

trim conditions used in this study are for conventional 

helicopters, the trimmed horizontal thruster position must be 
taken to equal zero. Minimum and maximum thrust perturbation 
movements are therefore the same as the total minimum and maximum 
movements . 


Finally, the user selects the auxiliary propulsion option by 
setting a flag either on or off. An off flag zeros the control 
derivative column in the B matrix (X^ — 0) . 

5.3.2 State Integration 

The state equation, x - Ax + Bu, is integrated forward in time 
from x using a fourth-order Runge-Kutta procedure. The 

O 

integration increment is constant and specified by the user. The 
stability and control derivative matrices are scheduled according 
to flight condition and thus are computed each time the integrator 
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moves to compute x. There are two entries in the state time 
history at interior point times. The second entry is necessary to 
account for discontinuities in x and 1 that occur at interior 
points (see Equation 5-43). 

5.3.3 Cost Calculation 

The integral of the cost function is computed by rectangular 
integration. The interior and final point costs are calculated 
using simple matrix mathematics. 

The cost of control alone is computed and retained. It is used 
within the optimization algorithm as a measure for evaluating 
whether the cost function has been sufficiently minimized. 
Sufficiency is established when the total cost has been reduced to 
a level where the cost of control makes up most (90%) of the total 
cost. Recall that we want to determine the control input 
histories necessary to fly the trajectory, and we are therefore 
not interested in restricting control excursions (K is small). 

u 

This means that when the trajectory is optimal, undesirable state 
excursions will be zero, and the (minimum) total cost will be 
comprised only of the cost of control, and thus the reason for 
using control cost as a minimization sufficiency measure. 

5.3.4 Adjoint Integration 

The adjoint equation (Equation 5-40) is integrated backward in 
time from l f (Equation 5-39) using a fourth-order, constant 
increment, Runge-Kutta procedure. Equation 5-40 contains the A 
matrix, which is a function of flight condition and therefore 
computed each time the integrator moves to compute 1. The 
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integration routine places two entries at interior point tiroes to 
account for discontinuities in 1 at t^s (see Equation 5-43). 
l(t) necessarily has the same number of time entries as x(t) and 
u(t). 

5.3.5 Control Saturation Adjoint Vector 

Each element of the adjoint, or influence, vector, m(t) , can 

take on positive or zero values throughout the interval. Elements 

of m equal zero for an unsaturated control. An element becomes 

positive when its corresponding control is saturated. When the 

control is saturated, its element in m(t) is solved using Equation 

* 

5-41. Because a fraction of this equation, i.e. the gradient H , 

* 

is used to perturb the control, and because — 0 when the 

control is saturated, any improvement in the cost by perturbing 
the control would violate the constraint. When (and if) the 

element of m Aegpioa to turn negative, this signifies that the cost 
can be improved (reduced) by wxbaixi/xxitLnx}. the control, i.e. , by 
perturbing the control by a fraction of the now nonzero gradient. 

5.3.6 Control Gradient Calculation 

t 

Equation 5-41 is used to calculate the gradient, H , for each 
scalar control. (When the control element is unsaturated, m^ - 

f 

0, and the gradients 3f and Ji are the same). The control 
gradient can take on positive, negative, and zero values at each 
(discrete) point in its time history, depending on which direction 
(positive or negative) and how far the control is from its optimal 
value. See the discussion in Section 5.3.7 below for its 
application. 
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5.3.7 Control Perturbation 


When the cost function is minimized, the gradient for each 
control equals zero over the entire interval, satisfying equation 
5-41. However,' to minimize the cost function requires numerical 
techniques because x(t) and i(t) are integrated in opposite 
directions in time, posing a two point bondary value problem. The 
steepest-descent method for solving this problem perturbs the most 
recent control history estimate by a fraction of the gradient 
history. The gradient gives the direction (positive or negative) 
and the oe£a£u$e magnitudes across time of the amount by which 
each control element should be perturbed. However, because this 
is a two point boundary value problem, which means the gradient 
expression is not "correct" until it exactly equals zero (which 
occurs when the cost is minimized and x(t), u(t), and l(t) are 
optimal), the gradient indeed provides only relative magnitudes, 
requiring that it be scaled by an "appropriate" fraction across 
the interval. As well, it cannot be expected that the same 
fraction is "appropriate" for each control element. Each 
element's fraction does, of course, apply over the entire 
interval . 

The "appropriate" fractions for the control gradient elements 
are the ones which generate a new (perturbed) control history, 
u (t) - u (t) - kH (t), which minimizes the cost within 

new previous u 

the pne&ent iteration of the trajectory optimization algorithm. 
[2,32] The search for the best fraction vector is an iterative 



procedure which takes place within each major iteration of the 
trajectory optimization algorithm. Its logistics for the three 
controls case are shown in the flow diagram below. The iteration 
takes place within the Perturb Control block of Figure 5-2. 
Component functions are explained following the diagram. 


Inputs from Main Program 



Perturb control ; integrate 
states; record new cost 








Figure 5-3 Flow Diagram for Fraction Vector Search 

5. 3. 7.1 Selecting fraction vector k 

To minimize the cost, we want to move the controls in the 
direction opposite the slope of their gradients [4] , by a fraction 
of their gradients; i.e., we want to subtract kH (t) . And since 

u 

the sign of the gradient tells us its direction, elements of k 
will always be positive (or zero) . 

Two options are available for selecting the fraction vector 
with which to begin this algorithm. One option is simply to use 
the vector which resulted from the previous search. A second 
method is employed when either: (1.) the current search is the 

first search within the first major iteration step (so that a 
fraction vector does not currently exist) or, (2.) the fraction 
vector from the previous search generates a control which yields a 
cost that is much greater than the current recorded cost. This 
second method takes, as the fraction for each control, 0.9% of the 





quotient of the average absolute value of the control over the 
entire interval, and the maximum of the absolute value of its 
gradient over the entire interval. 

5 . 3 . 7 . 2 Perturb Control 

Once the control is perturbed -- u (t) - u (t)- kH (t) 

new previous u 

-- the state integration and cost calculation are performed using 
exactly the same functions called by the main program. 

5. 3.7.3 Search Routine 

For m controls, the search routine searches the positive 
m- space by seaching along a starting search vector within this 
positive seach space, and then systematically searching along m-1 
orthogonal vectors . The search moves to an orthogonal vector once 
the cost-minimizing point along the current vector has been found. 
This procedure is common, and treated in many texts [2,32]. A 
simple example best illustrates how the routine works. 

Figure 5-4 plots cost contours versus fraction elements for a 
two control case. The fraction elements, k and which 

minimize the cost, are to be found. Starting at point A, vector 1 
is searched for the [k ,k ] pair which minimizes the cost along 
vector 1; this is achieved at point B. Then vector 2, which is 
orthogonal to vector 1, is searched, and the minimizing fraction 
pair is found at C. Then vector 3 is searched, et cetera, and the 
process continues until point [k is reached. 

To determine which direction to search along a given axis (from 
its nominal point) , the costs of points on each side of the 

nominal point, but very close to it, are computed. The direction 
of lesser cost is then searched. The routine defines two search 
lengths along every search vector, one on each side of the 
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starting point (see Figure 5-4), that extend from the starting 
point to the axes which the vector ultimately intersects, or to 
infinity. Costs are computed for four equally spaced points along 
this length. Of course, the search along the current vector stops 
when the cost increases, and if three points were investigated 
along that vector, a quadratic (least-squares) fit is performed to 
determine the point of minimum cost. An orthogonal search then 
proceeds from that point. 



Figure 5-4 Two-dimensional Gradient Search Vector Example 

Both the cost contours and the gradient change after each major 
iteration step, because x(t) and u(t) have changed. Therefore, 
the same point (set of fraction elements) may not reduce the cost 
on the next major iteration. Nevertheless, this point becomes the 
starting point for the orthogonal search. 

As one might expect, the greatest reduction in cost usually 
occurs after m vector searches, where m is the number of controls. 
Take, for example, Figure 5-4 again. In this two control case, 
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searches along the first two search vectors (vectors 1 and 2) 
reduce the cost more than all the subsequent searches. Because 
the integration and cost calculations take considerable time, and 
because they are performed for each point that is checked along a 
vector, the search routine checks to see if the cost after m 
vector searches is less than the starting cost (from the previous 
major step). If it is, the search ends. Though the cost may be 
further reduced by more searches, this cost reduction (usually 
small) is not worth the computation time, and it is therefore 
advantageous to recalculate the gradient (next major step) which 
should yield more substantial reductions in the cost. Of course, 
if the cost has not been reduced after m tries, another m vector 
searches are performed before going to the next major step. 

For this application, m - 3, and a vector search is carried out 
in the positive x-y-z space (or we could call it the k -k -k 3 
space) . The progression that this algorithm takes is to first 
search the r-z plane, which is the plane made by the z axis and 
the resultant vector r in the x-y plane . The second search vector 
is the vector orthogonal to the first vector in the r-z plane. 
There is an infinity of choices for the third vector, which is a 
vector orthogonal to the first two. The choice in this study is 
to use the vector in the plane parallel to the x-y axis at the 
point of intersection of the previous two vectors; i.e., at z - 
constant. 

Determining the * orientation of the intitial search vector is 
also a matter of choice. Eight distinct search vector 
orientations are possible in each plane, as shown in Figure 5-5. 
Unless a gradient element is zero over the entire interval, or a 
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control is not operating (e.g., when the auxiliary thruster is 
off), the vector from the origin (vector 1 in Figure 5-5) is taken 
as the initial search vector in each plane. In the special cases 
noted above, or in cases where the search along a vector 
terminates at an axis, the search vectors will be those vectors 
on, or parallel to, the representative axis. 

The algorithm codes the eight different vector types by their 
search lengths and angle of orientation at a given search point. 
The code changes systematically according to the orthogonalization 
of the vector. For problems with more than three controls, it may 
then become advantageous to replace this technique with a more 
general technique, such as a Graham- Schmidt orthogonalization 
method [ 2 ] . 

This procedure changes from one to three of the fractions each 
time the search routine moves to a new point. This tends to make 
the orthogonalization process and vector length computations 
somewhat complex. A parallel search system -- one that varies 
each fraction element independently, i.e., one that searches 
parallel to each fraction axis -- would be trivial, and greatly 
reduce complexity (and therefore some computation time) . This 
method was tried, but found to be inadequate. It is not as time 
efficient as the current algorithm. And, because it varies one 
element at a time, it tends to spend a lot of time in the "ridges" 


of the cost function contours. 




Chapter VI 

CONTROLLERS FOR THE AH-1G TRIMMED AT HOVER 

This chapter provides practical numerical examples of the 
development of the feedforward and feedback gains for an implicit 
and explicit model -following controller for longitudinal control 
of the AH- 1G in the hover flight regime. The ideal response model 
reflects the desired handling qualities criteria for a 
longitudinal velocity-command controller operating near hover in 
an NOE environment; numerical criteria data are taken from the 
references cited in Section 2.5. 

The development procedure includes the following analyses: 
input -output frequency response magnitude and phase analyses; 
eigenvalue placement and associated damping; and step input time 
history analyses. These analyses were used in designing each 
controller at each of the twenty-two trim points. Since it is 
rare that a controller can be developed which exactly meets the 
prescribed numerical handling qualities criteria, acceptable 
deviations, or "conditions of acceptability" were defined around 
these criteria. The controller was judged to be satisfactorily 
developed once the results of the frequency, eigenvector, and time 
history analyses met these conditions. These conditions also 
establish a common ground for comparing different controller 
structures. The conditions are defined in Section 6.4. Section 
6.5 discusses the methodology for selecting weighting matrices' 
elements and the sensitivity of the closed- loop dynamics to these 



elements. The examples below are for the hover trim condition, 
which was found to be the most difficult condition (of the 
twenty- two trim conditions investigated) for satisfying the 
handling qualities criteria. 


6.1 Dynamics of the AH-1G and the Handling Qualities Model 
The longitudinal model dynamics, x - Ax + Bu, are 
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It is obvious, in this simple model, that the u and w responses 
are completely uncoupled and have first-order responses with 
time - constants of — - — and seconds, respectively. q and 6 

.33 .40 

equal zero over all time since the perturbation initial condition 

vector, x , is a null vector and since q and 6 are unforced; this 
o 

means Q and 0 remain at their trim values, Q and 6 . 

o o 

The linearized longitudinal dynamics, x ■» Ax + Bu, of the AH-1G 
near hover (U-*1.69 ft/sec, W--.02 ft/sec, Q — 0 deg/sec, 0~-.73 

o o o o 

deg) are 
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( 6 - 1 ) 


The B matrix of the dynamics of the conventional helicopter has 
been augmented with a third column which incorporates the effects 
of the auxiliary propulsion system. The "2" in this third column 
is the control derivative for the auxilary propulsor's 
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contribution to the X- force (see Section 6.1.1 below). If the 
third column of B is null, meaning that no auxilary propulsion is 
employed, elements of the feedforward and feedback gain matrices 
corresponding to the thrust control will equal zero. This is 
shown by example for each controller. 

6.1.1 Inclusion of the Auxiliary Thruster in the Dynamics 

As discussed in Section 2.4, the advantages incurred using 
auxiliary propulsion come in the ability to achieve higher rates 
of acceleration and deceleration. We therefore derive the values 
for the thruster elements in the dynamic equation from the 
desired acceleration level, 0.3g to .35g. 

In this simplified representation of the auxiliary thruster, 
the thruster contributes only to horizontal longitudinal body-axis 
acceleration/deceleraticn, i.e., to u (see Equation 6-1). It was 
assumed that there are no changes in any of the stability 
derivatives or any of the control derivatives for the conventional 
helicopter; and, the first element in the third column of the 
augmented B matrix is the only non-zero element in that column, 
representing the control derivative 

The general expression for the dynamics of the thruster is then 
u - X T (6-2) 

T 

or, X t - u/T (6-3) 

We express the travel of the thruster control in inches, and 
assign its range to be commensurate with the ranges of the 
collective and cyclic: T can travel ±5 inches. For full control 
deflection at ±.3g, this yields 
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. 3g ( 32 . 174 ft/secVg) 
. I.93 = 2 


2 

s • in 


5 inches 

The AH- 1G ' s actuator lag time - constants for both the cyclic and 
collective inputs are 0.08 seconds; 0.08 seconds is also used for 
the auxiliary thruster's actuator lag time -constant . 


6.2 Implicit Model-Following Controller Design 

6.2.1 Designs Including Auxiliary Propulsion 

The controller gains for the AH-1G with auxiliary propulsion 
installed are derived first. 

The Q and R matrices that result in a contoller which exhibits 
the desired handling qualities are 
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The resulting weighting matrices are 
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The steady-state Riccati matrix, P, and influence matrix S are 
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The resulting feedback and feedforward 
gains are 
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-.0605 

-.0028 


.07365 

- .02042 

2.41 

-15.01 


.01263 

-.1946 


The resulting closed- loop dynamics are x - (A-BC^x + (BC^d: 


u 


-.3993 .00029 

.904 

.1716 


u 


.00185 

.3927 


u 












com 

w 


-.000028 -.33 

. 00176 

. 00267 


w 


.33 

. 0 00098 


W 


B 






+ 




com 

q 


.0326 -.0003 

- . 7821 

- .2928 


q 


- .00002 

- . 00045 


L J 

6 m 


0 0 

1 

0 


9 m 


0 

0 



The open- loop eigenvalues (A matrix) and closed- loop 
eigenvalues (A-BC^ matrix) are 

open-loop .1205 ±.2645i closed-loop -.3902 ± .37261 

-.^11 ±.19271 - .4019 

-.33 

Two of the open- loop eigenvalues are unstable. All four 
closed- loop eigenvalues are stable. Of the closed-loop 

eigenvalues, two have been properly placed at -.4 and -.33 on the 
real axis; they have been driven to the roots of the model. The 
two other eigenvalues have complex parts (f - 0.52); the frequency 
plots (Figures 6-1 and 6-3) and step input time histories (Figures 
6-2 and 6-4) show that they do not comprise the dominant modes of 
the closed- loop dynamic response, though they are nevertheless 
sufficiently damped. Inspection of the Bode plots and step input 
histories in Figures 6-3 and 6-4 show that the vertical response 
to a vertical command input exactly tracks the response of the 



first-order model to the same command; the responses of the other 
states are insignificant. 

The u to u response of the controller (Figures 6-1 and 6-2) 

com 

also emulates the model response very well; the horizontal 
velocity response, u, to a step input in u com tracks the model 
response almost exactly, and no pitching motion results from this 
command. The ability of the helicopter to translate forward 
without a nose -down attitude is attributed to ^the auxiliary 
thruster; for the conventional helicopter, discussed below, the u 
to u response is accompanied by a nose-down attitude. 

com 


6.2.2 Designs Without Auxiliary Propulsion 

For the conventional helicopter, the gains are computed in 
exactly the same manner as illustrated above, with the following 
two changes in the inputs: the open- loop B matrix is not 
augmented with the auxiliary propulsion dynamics, and the 
weighting matrices which ultimately provide the desired handling 
qualities are: 


Q 


900 0 100 100 

0 150 0 0 

100 0 0 0 

100 0 0 0 

- 


R - 


200 0 
0 10 


The resulting controller gains are 


c - 

.00955 

.00417 

-.0302 -.0477 

c - 

1 

.02595 -.0036 

1 

.2923 

-.0122 

.3504 -25.35 

4* 

.003917 -.257 


0 

0 

0 0 


0 0 


(The third row of gains contains zeros because the auxiliary 
propulsion is "off”). 
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The closed-loop eigenvalues are -.088 ± 2 . 016 . Again, two of 

-.4002 

-.3305 

the poles have been properly placed at the poles of the model. 
The frequency and step input response plots are shown in Figures 
6-5 to 6-8. Figure 6-5 shows a large amplitude response of q (and 
6 ) at a frequency of 2 rad/sec to a u input, and Figure 6-6 shows 

c 

the presence of the lightly damped mode in the pitch response to 
u inputs; the pitch rate feedback gain to u is only 0.35 
compared with -21 in the augmented helicopter's feedback matrix. 

This lightly damped mode does not meet the criterion set forth 
in Table 2-1. However, as noted above, it is not always possible 
to satisfy each criterion. In this example, increasing the 
damping of this mode (by increasing the values of the Q[l,3], 
Q[l,4], Q[3,l], and Q[4,l] weights) lengthened the time -constant 
of the horizontal velocity response to a horizontal velocity step 
input command, as shown in Figure 6 -6a. Because this study used a 
velocity-command controller, the velocity response criterion was 
judged to be a more important than the damping criterion. Section 
6-3 lists the order and importance of satisfying the handling 
qualities criteria that was followed in this study during the 
controller design. 

Comparison of the responses in Figures 6-1 through 6-8 for the 
conventional and thruster -augmented helicopters with implicit 
model -following controllers shows that both controllers exhibit 
good vertical velocity command responses. For horizontal velocity 
command inputs, the conventionally equipped helicopter has 
underdamped pitching motion as it accelerates horizontally away 
from hover. The controller for the augmented helicopter met all 
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the response criteria, and therefore provided decoupled horizontal 
and vertical velocity command responses and insignificant pitch 
response . 

6.3 Explicit Model-Following Controller 

6.3.1 Designs Without Auxiliary Propulsion 

Because the previous section showed the Ricatti matrix and 
closed loop matrices only for the controller design which included 
auxiliary propulsion, for the sake of completeness in showing the 
design synthesis the order of the previous section is reversed, 
and the controller gains for the AH- 1G mlthoui auxiliary 
propulsion installed are derived first. 

The Q and R matrices that result in a contoller which exhibits 
the desired handling qualities are 

Q * 20 o o o and R - 100 o 

0 20 0 0 L° 1 

0 0 0 0 

0 0 0 0 

The steady-state Riccati matrix, P , and influence matrices ? 2 
and S are 

- 28.65 
. 75468 
35.77 
943.87 

S * -.0172 .43669 

.19412 -.00404 
-.0719 .4337 
.07277 - 11.224 


r 1 

<w - ^ 


- 5.529 .07940 0 0 

.05084 - 3.333 0 0 

- 19.651 .16556 0 0 

18.07 .39948 0 0 


- 5.9488 

-.1244 

20.127 

- .1244 

3.3042 

- .3286 

20.127 

-.3286 

196.53 

- 28.65 

.75468 

35.77 
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The resulting gains are 


C « .00806 
1 

4,4409 

0 


-.41785 

-.10767 

0 


-.0042 -.0187 

-5.831 -42.94 

0 0 


.00742 

.42157 0 0 

C - 

-.0245 -.00062 

-3.975 

.07601 0 0 

3 

-.0106 

.49551 

0 

0 0 0 


0 

0 


The closed-loop 


dynamics are x 


(A-BC )x 


(BC )x + (BC )d: 

2 m 3 


u 


-5.7828 

.01328 

8.875 

23.48 


u 




w 


-.1098 

-5.6717 

.2564 

.17334 


w 




q 


.7217 

- .0195 - 

1.177 

-6.974 


q 




6 _ 


0 

0 

1 

0 


e. 






5,1534 

1.101 

.00894 

- . 2564* 


u 


-.00727 

.6424 


4. 

.76506 

56.261 • 

-.0006 

-.0033 


w 

+ 

.31036 

- .0078 


¥ 

-.64577 

-.0075 

-.001 

.00002 


q 


.00163 - 

. 08147 



0 

0 

1 

0 


e 


0 

0 



L 



- 


L J 

m 




U 

com 

W 

com 


The augmented dynamic equation is, 


X 


A-BC -BC ' 

X 


BC 1 


*■= 

1 2 


+ 

3 

X 

L m- 


0 F 

X 

L nr 


G 


and the closed- loop eigenvalues 
augmented state matrix f 


are 

A-BC 

l 

0 


the eigenvalues 


-BC 


The 


of the 
open- loop 


eigenvalues (A matrix) and closed- loop eigenvalues are 


open-loop .1205 ±.2645i 

-.4411 ±.19271 


Closed-loop -.2911 ± 1.8921 
-6.3726 
-5.6767 
-.40 
-.33 
0 
0 


The last four closed- loop roots are those of the prefilter (i.e., 
the model), and the first four are those of the closed-loop 
helicopter, all of which are stable. Bode plots are shown in 



Figures 6-9 and 6-11. The oscillatory modes (C - 0.15) can be 
seen in the time histories, Figures 6-10 and 6-12. Like the 
implicit controller for the conventional helicopter, this 
controller does not meet the damping criterion but was chosen over 
a controller which traded off greater damping for a longer 
horizontal velocity response time-constant (the time-constant was 
judged to be the more important of these two criteria for 
velocity -command controllers). Figures 6-11 and 6-12 show that 
the w to w response tracks the model response exactly. 

com 

6.3.2 Designs Including Auxiliary Propulsion 

X now has a non-zero value; the thruster is "on.” The Q 

T 

and R matrices are 

Q~ 10000 0 0 0 R * 2000 0 0 

0 500 0 0 0 2000 1800 

0 0 0 0 0 1800 2000 

0 0 0 0 

The resulting controller gains are 


C 

i 


-.1160 - . 4697 .1809 .4863 

-.8200 -.0016 - 22.59 - 36.47 

2.929 -.0273 15.45 14.28 


.11944 

.47337 

0 

0 

C * f- . 0245 
2 

-.00597 

1.1376 

.01843 

0 

0 

-.0082 

.350 

- 3.024 

.0097 

0 

0 

.0034 

-.109 


The closed- loop eigenvalues are - 2.03 ± 1 . 27 *. The dominant 

- 4.607 

- 6.4 

-.40 

-.33 

0 

0 

mode is complex and satisfies the damping criteria (f - 0.85). 
The other two roots have been driven far to the left along the 
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real axis. This controller satisfies the handling qualities 
criteria very well. The Bode plots and step input time histories 
for this controller are shown in Figures 6-13 through 6-16. The 
ability of the auxiliary thruster to reduce the forward pitch and 
oscillation of the helicopter out of hover is evident in the 
pitch rate/attitude time history of Figure 6-14. 

The explicit controller for the helicopter equipped with 
auxiliary propulsion provided u and w responses which track the 
model responses very well. However, though the controller met all 
the handling qualities criteria, it was not possible to design an 
explicit controller which equaled the ability of the implicit 
controller (for the augmented helicopter) to meet these criteria 
and also to provide unnoticeable pitching motion response to a 
horizontal velocity step input command. 

6.3 Conditions of Acceptability for Controller Analysis 

For the model dynamics of Section 6.1, the controllers of 
Sections 6.2 and 6.3 were judged to be acceptable and comparable 
in their abilities to satisfy the handling qualities criteria when 
they met the following conditions; these conditions were evaluated 
in the order given below: 

O The u/u and w/w frequency magnitude responses remained within 
c c 

+0.1 decibels of the model response from 0 to 1 rad/sec. 

O There were no peaks in these magnitude responses, and any 
valleys were confined to a 1 rad/sec span. 
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O The u/u and w/w frequency phase responses remained within ±4 

C C 

degrees of the model response from 0 to 1 rad/sec. 

O Any deviation in phase from the model was confined to -45 to 
-100° on the phase plot. 

O All input-output frequency magnitude responses other than u/u^ 
and w/w remained below -20 decibels for at least 5 rad/sec of the 

c 

0-6 rad/sec bandwidth and below -10 decibels over the entire band. 

O The time response of u to a step input in u^ and w to a step in 
w remained within il0% of the time responses of the model to 

c 

these same two inputs over the 5 second history. 

O The coupled responses, u to w and w to u , were confined within 
±0.05 ft/sec from zero over the 5 second time history. 

O Closed-loop eigenvalues were confined to the 0-6 rad/sec 
bandwidth . 

O The damping criterion (f > 0.5) was satisfied as long as it did 
not result in a system which violated any of conditions above. 
For both the implicit and explicit controllers without auxiliary 
propulsion, this criterion could not be met. This suggests that 
satisfactory damping using these controller structures is only 
possible when auxiliary propulsion is employed. 

O Pitch rate and pitch attitude time history responses were 
minimized while adhering to the conditions above. 

6.5 Selecting Controller Weighting Matrices. 

The elements of the Q and R matrices are parameters which must 
be chosen to yield a controller which provides the desired 
closed- loop properties. The following rules of thumb for 
selecting appropriate weighting matrices resulted from the 
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development of the controllers outlined in this study. Many of 
the rules discussed below are clearly set n in Figure 6-17, which 
plots root loci for the explicit controller (auxiliary propulsion 
installed) of Section 6-3 for various weighting matrices. 

O In general, the diagonal elements dictate the dominate response 
characteristics of the system. A few off-diagonal terms may 
occasionally be used to refine some aspects of the response, as 
discussed below. 

O The [1,1] and [2,2] elements of Q weight the u and w responses. 

They were varied independently to achieve the desired frequency 
magnitude and phase responses and time history responses. The 
magnitude of these elements varied from one controller to the next 
(from 20 to 10,000 in the examples in Sections 6.2 amd 6.3) and 
from one trim point to the next, but their influence in the 
frequency and time responses was easily discernible, so that the 
correct order of magnitude was quickly determined. These two 
weights affected only the response of their corresponding state 
variables, and once the correct order of magnitude was determined, 
changes within this order had virtually zero effect on the 
eigenvalue locations. These two parameters alone provided the 
desired decoupled velocity -command responses. 

O The [3,3] and [4,4] elements of Q showed little influence in the 
responses, except for the augmented helicopter's implicit 
controller. Occasionally, a unity weight was assigned to these 
elements to reduce pitching motions. For the implicit controller, 
at many trim conditions it was imperitive to weight the [1,3], 
[1,4], [3,1], and [4,1] elements of Q to reduce excessive 
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responses in q and 8 to u inputs. Once the [1,1] element was 
determined, these elements were found to be less than or equal to 
the [1,1] element. 

O The control weighting matrix, R, effectively controls the pole 
locations of the less dominant poles. In general, R did not 
strongly affect the dominant roots, and therefore it did not 
strongly influence the frequency and time responses of the system. 

In those cases where this influence was observed, the implicit 
controller was more susceptible to changes in the state responses 
from increasing R than was the explicit controller. 

O The unity matrix was a good choice for R. When (a) pole(s) were 
out of bounds of the 6 rad/sec frequency band, one or more of the 
diagonal elements of R was increased to pull the pole(s) back in 
bounds. A few tries showed which elements influenced which poles, 
and the magnitude of the proper element was adjusted to properly 
place the pole; often the poles were effectively placed by simply 
assigning the same scaling factor to all elements. 

O For the auxiliary thrusted helicopter, assigning a positive 
weight to the cross -coupled elements of R for the longitudinal 
cyclic and auxiliary thruster (elements [2,3] and [3,2]) provided 
an easy means of increasing the damping ratio of the dominant 
complex roots. The value of this weight was equal to or slightly 
less than, but of the same order of magnitude as, the smaller value 
of the [2,2] and [3,3] elements. 
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Figure 6-1 
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Implicit Controller: Hover, Aux Prop ON, u-command frequency response 
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Figure 6-2 
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Implicit Controller: Hover, Aux Prop ON, u-command step response 



Figure 6-3 
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Implicit Controller: Hover, Aux Prop ON, w-command frequency response 


w (model 



Implicit Controller: Hover, Aux Prop ON, w-commsnd step response 





Implicit Controller: Hover, Aux Prop OFF, u-command step response 
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Implicit Controller: Hover. Aux Prop OFF. u-commancJ step response 
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Figure 6-7 
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Figure 6-8 
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Implicit Controller: Hover. Aux Prop OFF. w-command step response 
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Figure 6-9 
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Explicit Controller: Hover. Aux Prop OFF, u-command frequency respon 
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Figure 6-10 
94 


Explicit Controller: Hover. Aux Prop OFF, u-corcmend step response 
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Figure 6-11 
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Explicit Controller: Hover. Aux Prop OFF. w-corcmand frequency response 
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Explicit Controller: Hover, Aux Prop OFF, w-command step response 
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Figure 6-13 
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Explicit Controller: Hover. Aux Prop ON. u-command frequency response 
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Figure 6-14 
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Explicit Controller: Hover, Aux Prop ON. u-command step response 
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Explicit Controller: Hover, Aux Prop ON, w-command frequency response 
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Figure 6-16 
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Explicit Controller: Hover, Aux Prop ON, w-command step response 
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Figure 6-17: Root Loci versus Q and R for Explicit Controller 

(auxiliary propulsion installed) 

Numbers refer to the pairs of Q and R matrices. 

' x ” denotes the helicopter's open-loop roots. 


Q and R matrices : 


1 


1 0 0 0 

0 500 0 0 

0 0 0 0 

0 0 0 0 


2000 0 0 
0 2000 0 
0 0 2000 


2 


500 000 

0 500 0 0 

0 0 0 0 
0 0 0 0 


2000 0 0 
0 2000 0 
0 0 2000 


1500 000 

0 500 0 0 

0 0 0 0 
0 0 0 0 


2000 0 0 
0 2000 0 
C 0 2000 


4 


3000 0 0 0 

0 500 0 0 

0 0 0 0 
0 0 0 0 


2000 0 0 
0 2000 0 
0 0 2000 


10000 0 0 0 


2000 0 0 

0 500 0 0 


0 2000 0 

0 0 0 0 


0 0 2000 

0 0 0 0 
»- - 



10000 0 0 0 
0 500 0 0 

0 0 0 0 
0 0 0 0 


5 0 0 
0 5 0 
0 0 5 


7 


10000 0 
0 500 

0 0 
0 0 


0 

0 

0 

0 


0 

0 

0 

0 


50 0 0 

0 50 0 

0 0 50 
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10000 000 
0 500 0 0 
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2000 0 
0 2000 
0 1000 


10000 000 
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0 0 0 0 
0 0 0 0 


2000 0 0 
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Chapter VII 

SELECTED NOE TRAJECTORIES 


Three representative longitudinal NOE maneuvers are defined 
in this chapter. The optimal control input histories for each 
maneuver are determined using the trajectory optimization 
algorithm. These trajectories, or maneuvers, are used in Chapter 
VIII to evaluate the controllers developed in Chapter VI, and to 
evaluate the use of longitudinal auxiliary propulsion. 

The collective and longitudinal cyclic ranges of stick travel 
for the AH-1G are [1]: 

6 : 0 to 10.7 inches B : -7.5 to 5.7 inches 

c Is 

7 . 1 Pop-up/Dash/Descent 

This maneuver, shown in Figure 7-1 is employed primarily when 
moving from one hiding point to the next. For example, the 
helicopter, hidden by tree cover, ascends above tree height, 
dashes across the tree -tops, and descends back to hover in a new 
clearing. The maneuver is done in minimal time. 
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Figure 7-1 

With auxiliary propulsion installed, the minimum time in which 
the maneuver can be effectively completed is 4.6 seconds. Table 
7-1 shows the specified values of the states at the interior and 
final point times, the trim states, the available control power at 
this trim condition, and the cost function weighting matrices. 
The specified control values (of the controls in the augmented 
state vector) at interior and final points are arbitrary since the 
corresponding weighting elements are zero; the weights nuuit be 
zero since we do not want to penalyze control excursions. 
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6 and B refer to stick trim; corresponding swashplate trim 
c Is 

o o angles are 14.83 and -.7 6 degrees. 

* ★ 

£ AY is defined positive upward 

Table 7-1 
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Figures 7-3, 7-4, and 7-5 show the spatial position and 
position histories, velocities, and control inputs for a 
conventional helicopter performing the pop-up maneuver. Figures 
7-6, 7-7, and 7-8 show the same plots for the auxiliary propulsion 
equipped helicopter. 

Auxiliary propulsion is a major advantage in this typ® of 
maneuver. The conventionally equipped helicopter has large 
nose-down attitudes; after one second of vertical climb, the 
helicopter pitches nose down to accelerate; the attitude returns 
to trim during deceleration (Figure 7-4). The path that the 
helicopter takes during the second half of its trajectory is not 
as square as it is for the augmented helicopter. "Cutting the 
corner" off this path is necessary to avoid excessively high pitch 
attitudes. If a more square path is desired, the forward speed of 
the helicopter must be further reduced before descending; this 
will, of course, increase the total time to complete the maneuver. 

This maneuver is a high agility task. As one expects, we find 
that the stick motions span their entire ranges (Figures 7-5 and 
7-8). Comparing these two figures shows the impact of the 
auxiliary thruster. For the conventional helicopter, the 
longitudinal cyclic saturates at its maximum during the dash, and 
then is pulled back strongly for the deceleration. For the 
auxiliary thrusted helicopter, the thruster saturates for 
acceleration and deceleration, and the amount of cyclic is 
reduced. The collective saturates for the thrusted helicopter 
because the maneuver is performed more quickly. 

The pitch attitude for the conventional helicopter is more 
extreme than for the augmented helicopter. The pitch rates for 
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both cases are understandably high, but fall within the acceptable 
limit of 0.69 rad/sec 2 : 0.56 for the conventional helicopter, and 

0.49 for the augmented helicopter. 

7.2 Bob-up at 40 knots 

The bob-up entails "hurdling" an object while flying at 
a moderate to high speed, as shown in Figure 7-2. 



Table 7-2 specifies the geometry and flying qualities for this 
maneuver. Both the conventional and auxiliary-equipped missions 


are accomplished in 2 . 5 seconds. 



interior and final points 


[ eux prop 
of f , on 


1.0 , 1.0 
1.5 ,1.5 


2.5 ,2.5 


trim conditions 


u w q 0 0 c Bis t AX AY 

ft/s ft/s deg/s deg/s in in in ft ft 


0 20 


ft/s ft/s 

67. 48 -2.13 


0 B 


3.75 -.0* 


X Y 

o o 

ft ft 


cost function weighting matrices 

State: diag [iso o 400 400 o o o iso o] 

control: diag [l l i] 

state -control: null matrix 

interior point weighting matrices 

diag [ 150 0 400 400 0 0 0 150 150 ] 


final point weighting matrix 

diag [ 150 150 400 400 0 0 0 150 150 1 


control power (inches) 

0 : - 3 . 7 , +10.0 B 


1# : -7 . 5, +5. 7 T: -5, +5 


e c and refer t0 8 tick trim; corresponding swashplate trim 
° ° engles are 12.77 and -.10 degrees. 


AY is defined positive upward 


Table 7-2 


Figures 7-9, 7-10, and 7-11 show the spatial position and 
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position histories, velocities, and control inputs for a 
conventional helicopter performing the bob-up maneuver. Figures 
7-12, 7-13, and 7-14 show the same plots for the auxiliary 
propulsion-equipped helicopter. 

Both the conventional and augmented helicopter were constrained 
to remain very near their trimmed horizontal velocity; this was 
effectively accomplished. The use of auxiliary propulsion 

decreases the amount of cyclic required; therefore, the pitch 
attitude for the augmented helicopter at the end of the maneuver 
is slightly less than the attitude for the conventional helicopter 
(the attitudes are well within the acceptable range). In either 
case, the pilot would have to reduce the pitch attitude at the end 
of the maneuver (alternatively, we could have weighted the 
constraint on the final pitch attitude more heavily) . There is no 
time advantage in using auxiliary propulsion here. There is a 
slight advantage in reducing the pitching motion during the 
maneuver . 

In both cases, the vertical velocity, W, is almost perfectly 
symmetric about its trim level over the interval, that is, it 
increases smoothly then decreases smoothly back to trim. The 
collective is increased, then decreased below its trim value to 
retard the vertical ascent, then increased above trim for the 
final half second to slow the descent speed to zero. 

2 

The maximum pitch rate never exceeds the 0.69 rad/sec limit; 

2 

the maximum rates (absolute value) are 0.27 and 0.4 rad/sec for 
the conventional and augmented helicopter, respectively. 
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7.3 Glideslope 

The glideslope is typically associated with descending from a 
moderate to high-speed cruise at an altitude of several hundred 
feet to a near-earth hover. The glideslope to ship deck landing 
is a common trajectory employing the glideslope. This study uses 
a very simplified glideslope. No particular glideslope path or 
angle is specified; the maneuver simply requires transition from 
flight at 40 knots and 250 feet altitude to a near-earth hover 250 
feet along the horizontal earth axis. 

Table 7-3 specifies the geometry and flying qualities for this 
maneuver. Figures 7-15, 7-16, and 7-17 show the spacial position 
and position histories, velocities, and control inputs for a 
conventional helicopter performing the glideslope maneuver. 
Figures 7-18, 7-19, and 7-20 show the same plots for the auxiliary 
propulsion equipped helicopter. 

The advantage of using auxiliary propulsion is evident in 
trying to fly the glideslope in minimum time while not exceeding 
the pitch attitude limits. The maneuver takes one second longer 
for the conventional helicopter than for the augmented helicopter 
(7 versus 6 seconds), and the pitch attitude is worse. In fact, 
as can be seen in Figure 7-16, the nose -up pitch attitude for the 
conventional helicopter at the end of the glideslope is not 
favorable; it exceeds the 10° limit. The spacial histories are 
nearly identical, as are the translational velocities histories, 
so the thruster provides the ability to rapidly decelerate while 
maintaining satisfactory attitude. 

The maximum absolute pitch rates over the trajectory are well 
within the allowable range, as should be expected. These values 
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7.4 Example of Path Evolution using the Trajectory Optimization 
Figure 7-21 shows three "paths" during different stages of the 
trajectory optimization for the pop-up/dash/descent maneuver (with 
auxiliary propulsion) . The initial control history yields a path 
that is very far from the desired trajectory. As the cost is 
minimized, the geometric specifications are met, as are the flying 
qualities criteria. The optimization of the other trajectories of 
this section converge similarly to their optimal path. 
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Chapter VIII 

EVALUATION OF IMPLICIT AND EXPLICIT CONTROLLERS 


The three maneuvers used for evaluation in this chapter are 
those maneuvers which were defined in Chapter VII. The optimal 
control input histories and state histories for these maneuvers 
(Chapter VII) are referred to often in the following discussion of 
the pilot command input histories. 

8.1 Pop-up Maneuver 

The pilot input histories for the implicit and explicit 
controllers are shown in Figures 8-1 and 8-2, respectively, for 
the conventional helicopter, and in Figures 8-3 and 8-4 for the 
helicopter equipped with auxiliary propulsion. A summary of the 
pilot workloads for this task is contained in Table 8-1. 

8.1.1 Without Auxiliary Propulsion 

For the implicit controller, the command inputs are choppy, and 
definitely coupled. Both stick motions reverse direction several 
times throughout the trajectory, and the rates of change are 
appreciable. Therefore, the pilot workload is very high. The 
explicit controller's command inputs are smooth. The vertical 
velocity command rises and then falls back to trim, in accordance 
with the desired shape of the pop-up and descent. The horizontal 
velocity command increases for the dash portion of the task, 
though it seems to remain unusually high at the end of the 
trajectory. 


8.1.2 With Auxiliary Propulsion 

The implicit controller, as for the case above, should be 
judged unacceptable. The explicit controller again accepts 
smooth, low pilot workload inputs. The horizontal velocity 
command this time, however, seems unusually low over the entire 
trajectory (Figure 8-4), especially considering the amount of 
positive cyclic and thruster input used throughout the maneuver 
(Figure 7-8). 

8.2 Bob -up Maneuver 

The pilot input histories for the implicit and explicit 
controllers are shown in Figures 8-5 and 8-6, respectively, for 
the conventional helicopter, and in Figures 8-7 and 8-8 for the 
helicopter equipped with auxiliary propulsion. A summary of the 
pilot workloads for this task is contained in Table 8-1. 

8.2.1 Without Auxiliary Propulsion 

The implicit model - following controller requires pilot inputs 
(Figure 8-5) that are not distinctively different from the 
swashplate inputs (Figure 7-11); there is a lot of movement in 
each stick, as witnessed in the command inputs and the rate of 
command inputs. Consequently, the pilot workload is relatively 
high. Furthermore, we notice that this controller has not 
decoupled the translational velocities very well. 

The explicit model -following controller, on the other hand, 
shows excellent controller characteristics. Figure 8-6 plots the 
pilot inputs using this controller. The vertical velocity input 
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is a very smooth bell -shaped curve. This command input reflects 
the curvature of the trajectory and the vertical velocity history 
(Figures 7-9 and 7-10); that is, the vertical response to the 
vertical command input is direct and predominately uncoupled. At 
the completion of the maneuver, when the helicopter is at its trim 
altitude and cruise speed, the W stick has been returned to its 

com 

trimmed location. Only minor horizontal velocity compensatory 
command is necessary, and this input is also very smooth. 

These smooth pilot inputs result in very low pilot workloads 
(see Table 8-1). The ability of the explicit controller to accept 
smooth inputs, mix them, and output the necessary "choppy" 
swashplate and thruster histories is quite pleasing. 

8.2.2 With Auxiliary Propulsion 

There is very little change in the pilot input histories for 
the thruster- equipped helicopter (because very little thruster 
power is needed to meet the specifications of this maneuver) . The 
same results concluded above apply to this controller. 

8.3 Glideslope Maneuver 

The pilot input histories for the implicit and explicit 
controllers are shown in Figures 8-9 and 8-10, respectively, for 
the conventional helicopter, and in Figures 8-11 and 8-12 for the 
helicopter equipped with auxiliary propulsion. A summary of the 
pilot workloads for this task is contained in Table 8-1. 

8.3.1 Without Auxiliary Propulsion 

The swashplate inputs for this maneuver are somewhat choppy 
(Figure 7-11). The pilot inputs for the implicit controller are 
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also choppy, whereas the explicit controller accepts smooth 
inputs. The smoothness of these inputs for the explicit 
controller, and the relative simple nature of this task results in 
very low pilot workload. 

8.3.2 With Auxiliary Propulsion 

Again, the implicit controller does not provide good controller 
characterisitics . The explicit controller is very ,£ood; the pilot 
histories differ from those above only in that the vertical 
velocity command at the end of the trajectory becomes positive to 
curtail the high descent rate which was made possible via the 
auxiliary thruster. 

In all of the cases presented above, the pilot workload for the 
explicit controller is considerably less than the workload for the 
implicit controlle The difference in workloads between 
maneuvers is attributed to the relative degree of agility inherent 
in the maneuver. 

Both controllers exhibited good frequency responses. It 
appears, therefore, that the explicit controller was able to 
accept smoother pilot inputs because of its use of a prefilter. 
The explicit controller attempted to match the actual and model 
state histories, whereas the implicit controller attempted to 
match the state rate histories; state rate matching required that 
the pilot input histories were similar to swashplate inputs so 
that the integrated state rates resulted in the correct state 
histories. Since the explicit controller matched state histories, 
the pilot inputs were required to satisfy the integrated model 


dynamics, and resulted in commands that were smooth and 
proportional to the desired velocities. 


Maneuver/ 
Aux Prop , 
Controller 


Pop-up/Off 

Implicit 

Explicit 

Pop-up/On 

Implicit 

Explicit 


Bob-up/Off 

Implicit 

Explicit 

Bob-up/On 

Implicit 

Explicit 


Glideslope/Of f 
Implicit 
Explicit 
Glideslope/On 
Implicit 
Explicit 


Pilot 

Workload 

Cost 


Open- loop 
Control 
Cost 


2 44 

8 66 

.0 32 


238 190 

238 190 




176 135 NA 

176 135 NA 

162 108 133 


21 162 108 133 


Normalized 

Cost 

W U 

com com 

e 

c 

B 

1 s 

.07 

1.05 

.02 

.20 

NA 

NA 

NA 

NA 

.09 

.74 

.05 

.23 

NA 

NA 

NA 

NA 

.06 

.14 

.01 

.23 

NA 

NA 

NA 

NA 
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Chapter IX 

CONCLUSIONS AND RECOMMENDATIONS 


Implicit and explicit model -following controllers have been 
investigated for use with the longitudinal degrees of freedom of 
an attack- type helicopter operating in a NOE environment. A 
computational trajectory optimization method was employed to 
evaluate these controllers, and to access the advantages incurred 
using auxiliary propulsion on helicopters operating in this 
environment. The auxiliary control power was sized according to 
desired acceleration potential. The contollers and longitudinal 
propulsion advantages were judged in terms of their resulting 
handling qualities, flying qualities, and pilot workload. 

A velocity-command controller, with decoupled command inputs, 
is preferred for these types of NOE tasks, which typically require 
high levels of helicopter agility. The explicit model -following 
controller uses a prefilter of model dynamics in the feedforward 
loop of the controller. This controller was superior to the 
implicit model -following controller in all cases tested, and was 
much easier to design because the same Q and R matrices were 
applicable for all trim conditions. Both could be designed to 
satisfactorily meet the prescribed handling qualities criteria 
using primarily diagonal weighting matrices. The explicit 
controller, however, required far less pilot workload. The 
explicit controller's command inputs reflected the direction, 
duration, and relative magnitude of the desired horizontal and 
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vertical velocity responses. This was true for the conventionally 
equipped helicopter, and for the thruster augmented helicopter, 
where the two command inputs were effectively mixed to provide 
inputs to the three helicopter controls (collective, cyclic, and 
thruster) . 

The auxiliary propulsor is very effective for NOE operations 
because (1) it increases the acceleration/deceleration potential 
of the helicopter, thereby decreasing the time to complete the 
task, and (2) the thruster- generated acceleration does not require 
excessive pitch motion, which is detrimental for navigation, 
target acquisition, and hiding. 

The trajectory optimization method provides a straightforward 
and efficient way to determine the optimal control input 
histories, which can be employed to then evaluate attributes of 
various controllers, propulsion systems, and perhaps other 
components of the helicopter. This study uses a soft constraint 
on the helicopter's pitch attitude and pitch rate. In retrospect, 
it may be more effective to use a hard constraint so as to 
absolutely limit the rate and attitude. In the current method, 
early stages of the optimization may yield rates and attitudes 
which are definitely unacceptable, and extreme attitudes may lead 
to divergence of the optimization. 

The controller gains are scheduled throughout the flight 
envelope. This is necessary to account for changes in the flight 
condition throughtout the trajectory. However, some gains change 
very little with flight condition, and some have negligible 
influence. This, in fact, was verified during early stages of 
this study, when the gains at the trimmed condition were used 


150 


throughout the entire trajectory. Therefore, further work in this 
area should be directed to correlating the most influencial gains 
with flight condition, and scheduling accordingly. For 
implementation in the aircraft, this will be necessary to reduce 
computational storage and real-time processing, and expense. 



Appendix A 

Helicopter Equations of Motion 

The most basic starting point in developing the dynamic 
equations of motion for the helicopter (or for any body) is to 
apply Newton's Second Law to the forces and moments acting on the 
helicopter, as described in reference 23. The helicopter is 
assumed to be a perfectly rigid body. Using a body-axis reference 
frame (Figure A-l), this yields the six degree of freedom 
expressions for helicopter motion. These equations relate the 
translational and angular accelerations to the aerodynamic, 
inertial, and gravitational terms. This study is concerned with 
the longitudinal equations of motion, and therefore the 
lateral-directional motion equations are not presented. 

Other assumptions employed in this derivation are that 

1. the helicopter has constant mass and mass distribution. 

2. the helicopter is flying at speeds low enough and 
distances short enough that the Earth can be considered 
flat; the body- axis frame then need not be related to 
rotation about the Earth's geocenter, and the inertial 
reference frame is therefore the flat Earth axis. 

3. quasisteady flow applies. 

The dynamics of the helicopter are presented in a body-axis 

frame. This is an axis system affixed to the body (Figure A-l). 

The X axis is forward through the nose of the helicopter; Y 
B B 

points out the right side; and Z fi points down through the body of 
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the helicopter. P, Q, and R designate rotational rates about the 
X , Y , and Z axes, respectively. The origin of this axis system 
is the helicopter's center of gravity. The XZ plane is assumed to 
be a plane of symmetry. Figure A-l also shows the 

Earth-referenced frame: points right; Y points out of the 

page; and Z^ points in the gravity (g) direction. Typically, it 
is necessary to consider the inertial to body transformation 
matrix which relates these two reference frames. -However, since 
we are dealing only with the linearized longitudinal motion 
equations, it is more succinct to state the simple result of the 
transformation linearization than to develop the transformation 
which is, for the most part, unused. (The result, q«0 , is 

included and explained below with the linearized longitudinal 
equations of motion) . 

The nonlinear longitudinal dynamic equations of motion are 

X - m(U + QW - RV + gsin0) - ma^ (A- la) 

2 - m(W + PV - QU - gcos©sin0) - ma (A-lb) 

z 

M - QI + PR (I -I ) - R 2 I + P 2 1 (A-lc) 

y x Z xz x 2 

m is the mass of the helicopter. Subscripted I and a are the 

moment of inertia and acceleration, respectively, along the axis 
(or in the plane) of the subscript variable. The equations 
represent the balance between the inertial and gravitational 
terms, and all aerodynamic and propulsive terms represented by the 
generalized symbols 2?, Z, and W. 

Since only the longitudinal degrees of freedom are of interest, 
these equations are rewritten without their lateral-directional 
terms (the lateral-directional terms are set to zero): 



2? - m(U + QW + gsin©) 


(A- 2a) 


Z - m(W - QU - geos©) (A-2b) 

W - QI (A- 2c) 

y 

These nonlinear equations are linearized by accounting for 
dynamic motion as small perturbations about an operating, or trim, 
point. The subscript ( ) denotes trim. We write the total 

o 

motion of each state as the sum of its trim value plus its 
perturbation value, i.e., we take the first two terms in a Taylor 
Series expansion. 
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In trim, nominal translational and rotation? 1 accelerations are 
zero: U - V? = Q - 0. Therefore, 

o O o 

u-u + u- u w-w + w- w Q-Q + q -= q (a-4) 

o o o 


and the trim equations are 
X - m(Q W + gsin# ) 

o o o o 

Z — m(Q U - gcos0 ) 

0 o o o 

M - 0 


(A- 5a) 
(A-5b) 
(A- 5c) 


Rewriting Equations A- 2 using Equations A- 3 and the derivative of 
Equations A-4, 

X + X - m[u + (Q +q) (W +w) + gsin(0 +0) ] 

o o o o 


- m[u + Q W +Qw+Wq+qw+ gsin# cos# + geos# sin#] 

o o o o o o 

(A-6a) 
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M + M « ql 

o y 


mfw • Q U • Q u • U q • qu - gcosfl cos0 + gsin0 sin0] 
0 0 0 0 o o 

(A- 6b) 
(A- 6c) 


This linearization technique assumes small perturbations, so 
that a small angle approximation is used: cos0 = 1 and sin0 £ 6. 

We also assume the products of perturbation variables are 
negligible. Using these two approximations, Equations A- 6 are 
X + X - m[u +QW +Qw+Wq+ gsintf + gcos0 6] (A-7a) 


Z + Z - m[w -QU -Qu-Uq- gcostf + gsin0 6] 

o o o o o o o 


M + M = ql 

o y 


(A- 7b) 
(A-7c) 


Subtracting the trim expressions, Equations A-5, from Equations 
A-7 yields the linearized (perturbation) equations: 

X - m[u + Q w + W q + geos# S] 


Z m[w-Qu-Uq+ gsin# 6] 

o o o 

M - ql 


(A-8a) 
(A- 8b) 
(A- 8c) 


or, 


Jx ■ 

-=u+Qw+Wq+ gcostf 6 

m o c> ° o 


(A-9a) 


2 

-«-w-Qu-Uq-f gsin0 6 (A-9b) 

m o o o 

Y - q (A- 9c) 

y 

Recall that X, Z, amd M represent the perturbation aerodynamic and 
propulsive forces and moment. 

The linearization of the inertial/body transformation yields 
the simple expression q«0 . This means that pitch attitude, 6 t is 
the direct integral of pitch rate, q. Thus, the independent 

perturbation variables in the longitudinal degrees of freedom are 



m[w -QU -Qu-Uq-qu- gcos0 co s9 + gsin0 sin 6] 

" 0 0 0 0 O o 


/A. 


M + M - ql (A- 6c) 

o y 

This linearization technique assumes small perturbations, so 
that a small angle approximation is used: cos0 = 1 and sin0 £ 8. 

We also assume the products of perturbation variables are 
negligible. Using these two approximations, Equations A- 6 are 



The linearization of the inertial/body transformation yields 
the simple expression q—£ . This means that pitch attitude, 8 t is 
the direct integral of pitch rate, q. Thus, the independent 
perturbation variables in the longitudinal degrees of freedom are 
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u, w, and q. With this in mind, we express X, Z, and M in terms 

of their sensitivities to changes in these independent state 

variables , and to changes in the longitudinal perturbation 

controls 6 and B (collective and cyclic, respectively): 
c Is 


3X , il ^ !X . ax 3X 

3^ u + + + art. 

c Is 


(A-lOa) 


_ az az az az . . az 

z “ au u aw w a q q ae e o aB B i s 

c Is 


(A-lOb) 


„ aM aM aM x 3H . x 5M R 

M “ aU u + a^ w + a^ q + JT e c + arts 

c Is 


(A-lOc) 


Dividing the force equations (X and Z) by m, and the moment 
equation (M) by I , the the resulting coefficients of the state 
variables are what are commonly referred to as the stability 
derivatives, and the control coefficents are the control 

derivatives. 

* _ A + I + i + A * + A dX E 

m in 3u m c)w m 9q m 86 c m 3B is 

^ c is 


tu + Xw + Xq + X- 6 + X B 

u w q pcc Bis Is 


(A- 11a) 


i az _i az i_ az i az. i_ sz 

m 3u U + m 3w W + m 3q^ m d$ c m 3B is 

C Is 


Zu + Zw+Zq + Z.0 + Z B 

u w q v c c Blsl 


(A-llb) 


M l 3M x l 3M i SH n l 3M . 1 3M 

I “ I a^I u I aw* I aq q I rt + I artis 

y y y y y c yls 


Mu + Mw + Mq + M. 6 + M B 

u v q u c c Bis Is 


(A-llc) 


Equating Equations A- 9 and A- 11, collecting like terms and 
rearranging, yields 


u - X u + (X -Q )w + (X -W )q - geos* 6 + X 6 + X B 

u wo qo o 0 c c Bis is 


(A-12a) 
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w ■* (Z +U )u + Z w + (Z +U ) q - gsin0 6 + Z 

o o w q o o 


q-Mu + Mw + Mq + M.0 + M B 

u w q B c c Bis Is 


* v + Z B 

pcc Bis Is 


(A- 12b) 
(A- 12c ) 


Writing Equations A- 12 in matrix form, and including the dynamics 
for 6 , we have a general expression for the linear longitudinal 
dynamics of the helicopter: 
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q 


X -Q 


Z +Q 
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u o 
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M. M 
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(A- 13 ) 


This matrix equation takes the form 

x - Ax + Bu (A- 14) 

where A is the stability derivative matrix, B is the control 
derivative matrix, x is the state vector, and u is the control 
vector . 
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q-Mu + Mw + Mq + M.0 +M B 

u w q^ 8c c Bis Is 

Writing Equations A- 12 in matrix form, and including 
for 6, we have a general expression for the linear 
dynamics of the helicopter: 



(A- 12b) 
(A- 12c ) 
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Scheduling 

The stability and control derivatives were scheduled according 
to the flight velocities U and W. Modal analyses [3,27] showed 
that the derivatives X , X , X , Z , X , and M , and 

u w q q Bis Bis 

variations in these derivatives, did not significantly affect the 
open- loop dynamic characteristics of the AH-1G. Thesefore, the 
average value of each derivative over the twenty- two trim 
conditions listed in Reference 1 was used throughout the entire 
flight envelope (see data above). 

The remaining stability and control derivatives were scheduled 
in a two-dimensional format by first curve - fitting the derivatives 
with respect to U, and then with respect to W, That is, they were 
scheduled as functions of U and W. Curve -fitting was either 
linear or quadratic, and the function for one derivative may have 
been composed of several distinct curves; all functions were 
continuous for all U and W. 

These non-constant derivatives were easily scheduled as 
functions of U by fitting curves to the data points. This method 
was also employed for scheduling with respect to W, except for the 
M , M , M , and Z derivatives. Because these derivatives varied 

u w q w 

considerably for different values of W, and because of the limited 
data set, curve-fitting was not directly possible. To facilitate 
scheduling, the effect of flight condition on these derivatives 
was rationalized by examining the equations that describe them in 
terms of their contributing aerodynamic forces and moments. By 
rationalizing the effects of flight velocities on rotor inflow, 
rotor angle of attack, thrust coefficient, and rotor flap, and 
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then looking at their effects on the derivatives, the derivatives 
were written as functions of U and VJ. 
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Appendix C 

Solution of the Steady-state Riccati Equation 


Reference [26] provides a detailed description of the following 
diagonal izat ion procedure for solving the steady- state Riccati 
equation. 

The equation has the form P — -PM -M T P + PNP + K and a 
solution is sought for P at steady- state , i.e., when P — 0. M, N, 
K, and P have dimension n x n. 

Define the 2n x 2n matrix 

7 [ M -N " 

* T 

K -M 


It is assumed that Z has 2n distinct eigenvalues. For each 
eigenvalue A, there is an eigenvalue -A; none of the eigenvalues 
have zero real parts. The eigenvalues of Z which are negative 
(stable) , are the eigenvalues of the steady-state closed-loop 
optimal controller. The diagonal modal matrix A contains the 
eigenvalues of Z with positive real parts. Z can be diagonalized 
as 



Partitioning W into n x n submatrices, W 


W W 
11 12 

w w 
21 22 


the steady- state solution of the Riccati equation is 


P 


w w ' 1 


22 12 


164 


ORIGINAL PAGE IS 
OF POOR QUALITY 

Appendix D 

Controller Gains at All Trim Points 

Trim numbers correspond to those in Appendix B. Numbers in the 
"Q" and "R" columns refer to the number which designates the Q and 
R weighting matrices used to develop the controller; these 
numbered matrices are listed at the end of the appendix. 


Implicit Controller -- conventional helicopter (no auxiliary prop) 
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Explicit Controller -- conventional helicopter (no auxiliary prop) 
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